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Abstract 

The extragalactic background light (EBL) is one of the fundamental ob- 
servational quantities in cosmology. All energy releases from resolved and 
unresolved extragalactic sources, and the light from any truly diffuse back- 
ground, excluding the cosmic microwave background (CMB), contribute to 
its intensity and spectral energy distribution. It therefore plays a crucial 
role in cosmological tests for the formation and evolution of stellar objects 
and galaxies, and for setting limits on exotic energy releases in the universe. 
The EBL also plays an important role in the propagation of very high en- 
ergy 7— rays which are attenuated en route to Earth by pair producing 7 — 7 
interactions with the EBL and CMB. The EBL affects the spectrum of the 
sources, predominantly blazars, in the ~ 10 GeV to 10 TeV energy regime. 
Knowledge of the EBL intensity and spectrum will allow the determination 
of the intrinsic blazar spectrum in a crucial energy regime that can be used 
to test particle acceleration mechanisms and VHE 7— ray production models. 
Conversely, knowledge of the intrinsic 7— ray spectrum and the detection of 
blazars at increasingly higher redshifts will set strong limits on the EBL and 
its evolution. This paper reviews the latest developments in the determina- 
tion of the EBL and its impact on the current understanding of the origin 
and production mechanisms of 7— rays in blazars, and on energy releases in 
the universe. The review concludes with a summary and future directions 
in Cherenkov Telescope Array techniques and in infrared ground-based and 
space observatories that will greatly improve our knowledge of the EBL and 
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the origin and production of very high energy 7— rays. 
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1. Introduction 



The extragalactic background light (EBL), defined here as the emission in 
the 0.1 to 1000 /im wavelength region, is one of the fundamental observational 
quantities in cosmology. It comprises the integrated light from resolved and 
unresolved extragalactic sources, and the light from any truly diffuse back- 
ground, excluding the cosmic microwave background (CMB). It is therefore 
the repository of all energy released by nuclear and gravitational processes 
since the epoch of recombination. A significant fraction of this radiation is 
shifted by cosmic expansion and by absorption and reradiation by dust into 
infrared (IR) wavelengths. Consequently, its intensity and spectral shape 
hold key information about the formation and evolution of galaxies and their 
stellar and interstellar contents throughout cosmic history. A strict lower 
limit on the EBL intensity is provided by the integrated light from resolved 
galaxies, hereafter referred to as the integrated galaxy light (IGL). 

The EBL plays also an important role in the propagation of high energy 
7— ray rays that are predominantly emitted by blazars, a subgroup of active 
galaxies hosting active galactic nuclei (AGN), whose relativistic jet is pointed 
towards the Earth. High energy photons emitted by blazars are attenuated 
by photon-photon interactions with the EBL, a process that can be used to 
set important limits on both, the intrinsic spectra of blazars and the intensity 
of the EBL in select energy and wavelength regions where these interactions 
are most prominent. 

The EBL is intimately connected to the diffuse X-ray, radio, and super- 
nova neutrino backgrounds. Deep X-ray surveys have resolved the X-ray 
background into point s ources, most of which are dust enshrouded AGNs 
( IMushotzky et al.l . |2000| ). Up to 90% of the X-ray energy produced in in- 
dividual AGN ca n be degraded and reradiated predominan t ly at mid-IR 



wavelengths (e.g. iFranceschini et al.l . |2002| ; iBallantyne et al.l . 120061 ) . Con- 



sequently, the X-ray background can be used to predict the EBL intensity at 
at these wavelengths. Current estimates s how that about 15% of the 24 /xm 



EBL intensity is powered by AGN activity ( jTreister et al.l . 120061 ; ISoifer et al. 
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20081 . and references therein). Conversely, the connection between mid-IR 



bright sources and AGN can be used to estimate the contribution of ob 



scure d AGN to the X-ray background (IGandhi fc Fabianl . 120031 ; ISoifer et al. 



2008 



and references therein). 
Massive stars that power the IR emission also emit radio free-free emis- 
sion during the main sequence phase, and radio synchrotron emission during 
the supernova remnant phase of their evolution. The IR emiss ion from star- 



forming galaxies is therefo re correlated with the radio emission ( iLisenfeld et al. 



19961 : ICondon et al.l . Il99ll ). This correlation can be used to estim ate the con- 



tribu t ion of star-forming galaxies to the cosmic rad io background (IHaarsma fc Partridge] . 



19981 : iDwek & Barked . 120021 : IPonente et al.l . 1201 if ). 



Most of the EBL intensity is powered by massive stars that end their 
life as core collapse supernovae. The total EBL intensity can therefore be 
used to derive an es timate of the supernova rate and th e resulting flux of 
supernova neutrinos (IHoriuchi et al.l . 120091 : iBeacoml . 120101 ) . The detectability 
of these neutrinos can be greatly enhanced by the proposed introduction 
of gadolinium in existing large water Cherenkov detectors (such as Super- 
Kamiokande) (IBeacom &: Vagina . 120041 ) . Gadolinium has a very high capture 
cross section for neutrons generated in v e + p — > e + + n reactions, and can be 
introduced in the form of soluble trichloride (GdCL>). Following the neutron 
capture, the Gd emits an 8 MeV 7— ray which produces relativistic electrons 
by Compton scattering. The Cherenkov radiation from these electrons is 
more easily detected than that produced in the cascade of the 2.2 MeV 
7— ray generated by the capture of neutrons by free protons. 

Several reviews have appeared in the literature, presenting a histori- 
cal overview of the importance of the EBL, early estimates of its inten- 
sity, the quests for its detection, and it s many astrophysical implications 



( IHauser &: Dwekl . l2001t iKashlinskyl . l2005t lLagache et al.l . |2005|). Since these 



reviews were written, significant advances have been made in studies of the 
EBL with the launch of UV (Galex) and IR space observatories (Spitzer, 
Herschel, and Akari). These observatories, together with ground-based tele- 
scopes, such as 2MASS, have provided new limits on the EBL ranging from 
UV to submillimeter wavelengths. Deeper galaxy number counts and new 
data analysis techniques of stacking astronomical images have narrowed the 
gap between the contribution of resolved galaxies and the true intensity of 
the EBL. 

The Fermi Gamma-ray Space Telescope, operating between 200 MeV and 
300 GeV, and ground-based air Cherenkov detectors (H.E.S.S., MAGIC, and 
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VERITAS) operating in the ~ 50 GeV to 100 TeV range have broadened 
the energy window for the studies of 7— ray sources. These advances have 
led to the detection of new GeV and TeV 7— ray sources and provided new 
data for dete rmining their i ntrin sic spectra. Reviews of the se subjects were 
presented bvlWeekesi (120081 ) and iHinton &: Hofmannl (120091 ) . More recently, 
Dermerl (I2OI2F presented a review of the Fermi catalog of 7— ray sources and 



the physics of the production of relativistic particles and 7— rays from these 
sources. Table 1 presents a glossary to the acronyms of the observatories and 
instruments referred to in this review. 

These developments provide the main impetus for this review. We first 
present, in §2, the basic formulae describing the attenuation of photons by 
pair producing interactions with other photons. We then show how this at- 
tenuation will affect 7 rays traversing a radiation field characterized first by 
a pure black body, representing the stellar emission component of the EBL, 
and then by a more realistic EBL that includes the dust emission component. 
This attenuation can, in principle, be used to determine the intensity of the 
attenuating radiation field if the intrinsic source spectrum is known. In §3 we 
survey the type of 7— ray sources that are used in these studies, their spec- 
tral characteristics, the physical mechanisms for generating their spectra, and 
constraints on their spectral shape imposed by general physical principles. In 
§4 we summarize measurements and limits on the EBL intensity determined 
by direct measurements and by adding the light from resolved galaxies. Mod- 
els for the EBL intensity and its evolution with redshift are summarized in 
§5. In §6 we summarize the constraints on the EBL intensity derived from 
7— ray observations of blazars, emphasizing the different assumptions made 
on the intrinsic blazar spectra to derive these limits. EBL models predict 
the 7— ray opacity of the universe at different energies, and in §7 we compare 
these model predictions with blazar observation. Throughout this review it 
was tacitly assumed that the production of 7— rays takes place exclusively in 
the sources. In §8 we consider alternative scenarios of 7— ray production that 
could have important implications for EBL limits, namely, that a significant 
fraction of the observed 7— rays could be produced en route to Earth. The 
role of the EBL in setting limits on exotic energy releases in the universe 
in briefly discussed in §9. A summary and future prospects for the fields of 
7— ray and EBL research is given in §10. 
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2. The EBL and the Attenuation of Gamma-Ray Photons 



2.1. The EBL 

The differential specific flux at wavelength A , dF u (\ Q ), received from 
radiative sources within a co moving volume e lement dV c (z) at redshift z at 
wavelength A is given by (e.g. 

dF„(A ) 



Mo et all 12010( 1: 



'l + z. 



£ u {\, z) dV c (z) 
4?r d L {z) 2 



(1) 



where C U (X, z) is the comoving specific luminosity density of the sources, di, 
is their luminosity distance, and the (1+z) factor arises from the decrease in 
energy of the emitted photons due to the redshift, and A = (1 + z)\. 

The specific comoving intensity of the EBL per unit solid angle, 5Q, at 
redshift Zq and wavelength Ao is given by an integral over all energy releases 
over cosmic history: 
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where c\dt/dz\ is given by (e.g. Mo et al.L 120101 ): 
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H is the Hubble constant, and Qr fl m , Qf~ and are the dimensionless 
density parameters of the radiation, matter, the curvature, and the cosmo- 
logical constant A, obeying the relation: Qr + Q m + + Q\ = 1. The second 
expression for E(z) is for a matter dominated (Qr << 1) universe, and the 
third is for one that is matter dominated and flat (Qj. = 0). In the concor- 
dance cosmology mod el: H = 70 km s -1 Mpc -1 ; Q m = 0.27, and f2 A = 0.73 
flHinshaw et all l2009h . 
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2.2. Gamma-ray attenuation by pair production 

The interaction between two photons with energies and e^, will lead 
to the creation of a particle anti-particle pair when the total 7— ray energy 
in the center of momentum of the system exceeds the rest frame energy of 
the two particles. The threshold for the creation of an e + +e _ pair is given 
by: 

2 (m e c 2 ) 2 

eth {R r ^z) = ^- ir) (5) 




Figure 1: Schematic illustration of the 7 — 7 pair production reaction, showing the defini- 
tion of the angle 8 between the interacting photons. 

where \i = cos^, and 6 is the angle between the two photons, as illustrated 
in Figure [U 

The cross-section for the 7 — 7 interaction is given by: 



cx 77 (£ 7 , e, fi, z) = (1 - f3 2 ) 
where 



2(3 {(3 2 - 2) + (3 - /3 4 ) In 



(1-/3) 



(6) 



^V^-t) (7) 

Figure [2] (left panel) depicts the cross section as a function of (3. The 
cross section peaks at a value of (3 = 0.70, providing a relation between the 
energies _E 7 and e (or wavelength A) at the peak, given by: 

The right panel of the figure depicts the cross section as a function of b = 
2(mc 2 ) 2 / E^e for different values of the angle 9. When the photons are moving 
in the same direction (6 = 0), the cross section collapses to a delta-function 
at b = 0, and the energy threshold becomes infinite. 



6 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.5 1.0 1.5 2.0 

/S b = 2(m e c z ) z /E,€ 



Figure 2: The cross section for the 7 — 7 interaction. Left panel: its dependence on j3 
[eq. (7)]; Right panel: its dependence on b for different angles of incidence. 

2.3. The Attenuation of 7— rays from Cosmological Sources 

En route to Earth, 7— rays from cosmological sources have to pass through 
the radiation field of the EBL, resulting in their attenuation by pair producing 
interactions. The optical depth of a 7— ray photon at an observed energy E y , 
emitted by a source at redshift z due to this process is given by: 

f z di f 1 1 — a f°° 
T 77 (£ 7 ,z)=y dz' — J dn^^- J den^z'^l + z'fa^'.z') (9) 

where n e (e, z) = dn(e, z)/de is the specific comoving number density (cm -3 eV _1 ) 
of background photons with energy e at redshift z, and the (1 + z) 3 term 
represents its conversion to a proper number density. The pair-production 
threshold energy is e' th = 2(m e c 2 ) 2 /E 1 (l — + z), where the (1 + z) fac- 
tor takes into account that the observed 7— ray photon had a higher energy 
at the redshift of the interaction. The parameter ft' — (1 — e^/e) 1 / 2 , and 
d£/dz = c\dt/dz\, where £ is the proper distance. 

Calculating the EBL opacity to 7— rays from cosmological distant sources 
requires knowledge of the evolution of the comoving specific photon number 
density n e (e, z) as a function of redshift. The specific number density of 
photons with energy e at redshift z is related to the specific EBL intensity 
at a given redshift z by: 

4:71 

e 2 n e (e,z) = — v I u (u, z) (10) 
c 

= 2.62 x 10~ 4 z/^(z/,2) 

where e = hv, I v (u,z) is given by eq. (2), and the coefficient in the second 
line was calculated for e in eV, n e in cm -3 eV -1 , and v I v in nW m~ 2 sr -1 . 
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Finally, we point out that the 7 — 7 cross section is wide, so that in 
calculating the 7— ray opacity, strong variations in the EBL spectrum are 
smoothed out over a wide range of 7— ray energies. The EBL intensity at a 
given wavelength is therefore effecting r 77 over a wide range of 7— ray energies 
around the peak given by eq. (8). 

2.4- A Simple Example: An EBL given by a diluted blackbody spectrum 

Of particular interest is the behavior of r 77 for a background radiation 
field that is represented by a diluted blackbody. Figure [3] (upper left panel) 
depicts a local EBL characterized by a Planck function, normalized to an 
intensity of 10 nW m~ 2 sr -1 at 1 /im. The upper right panel of the figure 
depicts the photon number density. The bottom left panel shows the 7— ray 
opacity at redshift z = 0.2, assuming a non-evolving EBL, and the right 
panel shows the source attenuation as a function of 7— ray energies. Also 
shown in the figure are the energy regimes in which substantial changes in 
the slope of the opacity occur (dashed lines). 

The rapid rise in the EBL spectrum between 0.5 and 1 /xm results in a rise 
of the 7— ray opacity, and the onset of substantial source attenuation in the 
10 to 500 GeV energy region. This sudden increase in the GeV attenuation 
creates a break, TceV, in the spectrum, defined as the difference in power law 
index between the unattenuated and the attenuated region of the spectrum 
(see Figure in this paper). At higher 7— ray energies, the spectrum of a 
blazar characterized by an intrinsic power law will exhibit a second spectral 
break around ~ 1 TeV. For an evolving EBL, the magnitude and location 
of this spectral break are expected to evolve with redshift. The substantial 
decrease in the attenuation at a few TeV is a consequence of the particular 
choice of the EBL spectrum, which decreases rapidly at wavelengths beyond 
~ 2 fim. 

2.5. A More Realistic Example: An EBL that includes dust emission 

Figure H] depicts a more realistic presentation of the current EBL spectrum 
(left panel) and the 7— ray opacity for d i fferen t redshifts (right panel), taken 



from model calculations of iFinke et al.l (120101 ). At wavelengths short wards 



of ~ 5 yum the spectrum represents the stellar and AGN contributions to the 
EBL. At longer wavelengths the spectrum represents the AGN and starlight 
energy that was absorbed and reradiated by the dust. The right panel shows 
the energy dependence of the 7— ray opacity for sources at different redshifts. 
The opacity calculations took into account the evolution of the EBL with 
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z = 




0.01 0.10 1.00 10.00 100.00 0.01 0.10 1.00 10.00 100.00 

E 7 (TeV) E T (TeV) 

Figure 3: Top left: A diluted black body representation of the stellar emission component 
of the EBL; Top right: The corresponding proper photon number density versus energy; 
Bottom left: The 7— ray opacity versus energy, E~; Bottom right: The 7— ray attenu- 
ation. The figure illustrates the dramatic change in the attenuation at the 7— ray energy 
that corresponds to the wavelength at which the slope of the EBL spectrum changes. The 
different slopes are depicted as dashed lines in the figure. 



redshift. The figure illustrates the relation between the EBL spectrum and 
the energy dependence of the 7— ray opacity. The initial rise of the EBL 
intensity at UV-optical wavelengths causes an increase in the 7 — 7 opacity 
between 10 and 500 GeV. The decline in the EBL intensity between ~ 1 and 
15 jum causes r 77 to rise less rapidly between 1 and 10 TeV. The slope of r 77 
in this region reflects the ratio of the ~ 1 to 15 /im intensities of the EBL. 
The rise in r 77 beyond 10 TeV reflects the rise in the EBL towards the peak 
of the dust emission at ~ 100 — 200 /im. 

The energy dependence of r 77 will give rise to several breaks in the spec- 
trum of 7— ray sources that reflect the changes in the slope of the opacity. 
The first spectral break, ATceV occurs between 10 and 500 GeV. The second, 
AT-TeV around 1 TeV, and the third around 10 TeV. 

The first break has been used in most EBL studies to date with various 
assumptions on the intrinsic source sp ectra, and t he se cond break has been 
most recently explored in the analysis of lOrr et al.l (120111 ). A review of studies 
that utilize the first and second break for constraining the EBL is presented 
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Figure 4: Left panel: Calculated EBL intensity versus wavelength at z = ; Right 
panel: The 7— ray opacity versus energy for sources at different redshift (see labels). The 
figure illustrates the correlation betwee n the change s in th e slope of the EBL intensity 



with those in r 77 . Model calculations by iFinke et al. I (|20ld) . Details in §2.5 of the text 



in §6. A break at ~ 10 TeV has yet to be discovered. Starburst galaxies, 
which have a hard 7— ray spectrum, are the most promising subject for such 
analysis. 

3. The Types and Spectra of Extragalactic GeV/TeV Sources 

Determination of the EBL intensity from GeV-TeV 7— ray observations 
requires knowledge of the intrinsic spectrum of the sources. Here we list 
the different sources, their spectral characteristics, the different proposed 
mechanisms for their 7— ray production, and the physical limits on their 
spectral energy distribution at very high energies. 

3.1. The Types of Extragalactic GeV/TeV Sources 

The currently available GeV-TeV 7— ray sources that are being used to 
derive limits and constraints on the EBL are listed in Table 2. They in- 
clude the accretion-powered relativistic jets of active galactic nuclei (AGNs), 
namely BL Lacertae objects (BL Lacs), flat spectrum radio quasars (FSRQs) 
and a few radio galaxies. The list is complemented by the recent detections 
of two nearby starburst (SB) galaxies. In contrast to AGNs, their 7— ray 
spectrum is generated by the cumulative effects of cosmic-ray acceleration 



in sh ocks generated by a large number of supernova remnants (IVolk et al 



19961 ) . 



Blazars : To date, the most numerous sources used in EBL studies are 
blazars. Historically, they have been divided into two sub-classes based on 
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their optical properties: FSRQs, characterized by strong emission lines; and 
BL Lacs, characterized by weak or lack of emission lines. Because of the weak- 
ness of their emission lines, the redshift determination of BL Lac blazars has 
proven difficult or even impossible in many cases. The status of blazars as 
bright GeV and TeV sources arises from the fact that their relativistic jets are 
closely aligned with the observers line of sight. Consequently, the luminosity 
of a 7— ray emission region moving relativistically along the jet axis in the 
direction of the observer is strongly beamed, enabling its detection at cos- 
mological distances. Occ asional strong flaring ac tivity rend ers the following 



BL L acs: PKS 2155-304 (lAharonian et al.l. l2007fi); Mrk 501 (ICatanese et al 



1997: 



Abd o et all. 12011 a!) and Mrk 421, flGaidos et all . Il996l : lAcciari et al 



2011al : lAleksic et all 1201 II) the brightest T eV sources; a nd the following FS 



RQs: 3C 454.3 flDonnarumma et all 120091 ). and 3C 279, flWehrle et allll998h 



the brightest GeV emitters in the sky. The flaring has provided high quality 
7— ray spectra and has led to their detection at redshifts as far as z ~ 0.5 at 
TeV energies with IACTs, and as far as z ~ 3.2 at ~ 10 GeV energies with 
Fermi. 

The combined GeV- TeV observations of blazars make it possible to study 
their spectra over a larger range of redshifts, thereby enabling the studies of 
the EBL over a wider range of wavelengths. GeV photons interact mainly 
with UV/optical photons, whereas TeV photons probe mainly the near- to 
mid-IR region of the EBL. Since the intensity of the EBL is much lower at 
UV energies, the universe is transparent to 7— rays below 10 GeV, becoming 
essentially opaque for TeV sources at redshifts of z > 0.5. The Fermi Gamma- 
Ray Space Telescope provides important probes of the UV region of the EBL, 
and the GeV transparency can be used to test evolutionary models of the 
EBL to relatively large redshifts (z > 1). 

Radio galaxies: The jets in radio galaxies are significantly misaligned with 
respect to the observer's viewing direction, and thereby provide no relativistic 
Doppler boosting. This limits the detection of radio galaxies with current 
generation 7— ray telescopes to the local group and the Perseus galaxy cluster. 
Deep 7— ray observations of radio galaxies with CTA combined with spatially 
resolved studies in the radio, optical and X-ray will play an important role in 
understanding the physics of relativistic jets. These observations are likely 
to yield spectra up to w 10s TeV which will provi de useful con s traints on 
the EBL in the mid- and far-IR wavelength regions (lAcciari et all 12009a . for 
the VERITAS Collaboration, the VLB A 43 GHz M 87 Monitoring Team, the 
H.E.S.S. Collaboration, and the MAGIC Collaboration). 
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With sufficiently high spatial resolution, the 7— rays produced by IC scat- 
tering of CMB and EBL photons off the relativistic electrons of the lobes of 
radio galaxies can be used to s et limits on their energy dens ity in the imme- 



diate vicinity of these objects (jGeorganopoulos et all 12008k se e §6 below) 



Starburst galaxie s : The detection of starb urst galaxies M82 (Acciari et al 



2009bl; lAbdo et all l2009h and NGC 253 flAcero et all |2009|), potentially 



opened a new wavelength regime for studying the EBL. The 7— rays in star- 
burst galaxies are generated by cosmic rays that are accelerated by a large 
number of supernova remnants, giving rise to hard 7— ray spectra that ex- 
tend to energies of 10s of TeV. The ~ 10 TeV opacity to nearby starburst 
galaxies is quite small, and about unity at energies of ~ 50 - 100 TeV. Nearby 
starbursts are therefore important probes of the EBL at far-IR (~ 100 /im) 
wavelengths that cannot be probed by other 7— ray sources because of the 
relative softness of their spectra compared to those of SB galaxies. 

3.2. The Spectra of Extrag atactic GeV/TeV Sources 

Over a sufficiently small energy range the blazar spectrum can be charac- 
terized by a power law, dN/dE oc E _r , with different indices, TceV and rVeV, 
at GeV and TeV energies, respectively. An important characteristic of the 
observed spectra is the presence of a break, defined as ALcev = TceV — TTeV, 
occurring between GeV and TeV energies, the exact location depending on 
the source's redshift. A source with an intrinsic spectrum characterized by 
a single power law out to TeV energies will have a value of ALgcv — 0. 
Without any intergalactic absorption this value will remain constant with 
redshift. 

The spectral index Tcev is obtained from a power law fit to the ~ 1 — 
10 GeV region of the spectrum which is unaffected by EBL absorption. If 
the intrinsic blazar spectrum is an extension of this power law to energies of 
~ 1 TeV, then any spectral break {Yrev > ^Gev) in the observed spectrum 
can be regarded as evidence for EBL absorption. 

A spectral break analysis of the amount of EBL absorption provides there- 
fore a powerful method for studying the EBL. It is a differential method that 
replaces knowledge of the intrinsic blazar spectrum with a weaker require- 
ment, namely that the power law representing the intrinsic blazar spectrum 
at GeV energies can be extended to TeV energies as well. 

Table 2 lists the values of Tcev an d ^TeV and the redshifts for all GeV and 
TeV detected blazars. Almost all sources exhibit a spectral break (Arcev < 
0) at energies between 10 GeV and 1 TeV. Figure 5 depicts the dependence 
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Figure 5: The difference between TQeV, the spectrai index at GeV (Fermi) energies, 
and r Tc v, the energy spectral index in the TeV regime (H.E.S.S, MAGIC, VERITAS) 
is shown as a function of their redshift. Red squares (radio galaxies), red stars (starburst 
galaxies), empty circles (HBLs, high-frequency peaked BL Lacs), blue downward triangles 
(intermediate- frequency peaked BL Lacs), filled circles (LBLs, low- frequency peaked BL 
Lacs), red upward triangles (FSRQs, flat spectrum radio quasars) indicate the different 
types of 7— ray sources. 

of Arcev on redshift. The figure shows a clear trend of increasing |Ar"G e v| 
with redshift, strongly suggesting that the break is the consequence of the at- 
tenuation of the source spectrum by the EBL. As the optical depth increases 
with redshift, the observed 7— ray spectrum becomes softer, the position of 
the break moves to lower energies, and Aro e v becomes more negative. The 
detailed redshift dependence of ArQev reflects the evolution of the spectrum 
and proper photon number density of the EBL with redshift (see §7). 

Figure also shows that there is significant scatter in ArQ e v & t an Y given 
redshift. This suggests that some sources have considerable intrinsic harden- 
ing/softening in their spectra. Indeed, observations show that sources exhibit 
a wide range of spectral trends, i.e., spectral softening (3C 279, PKS 1510- 
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08) an d spectral hardening (1ES 0502+675) in the GeV regime (lAbdo et al. . 



2010bl ). while the spectra of other blazars (Mrk 421, lAbdo et all l201ld) 



(Mrk 501. lAbdo et alllioilah extend without any cut off to 300 GeV flAbdo et al 



2010a), or multi-TeV energies ( IKrennrich et all 120021 ). depending on the flar- 
ing state of the source. A lower limit in the scatter of Arcev at a given 
redshift may suggest the combined effects of intrinsic spectral softening of 
the source spectrum and the effects of EBL absorption. An upper limit in 
the scatter may simply indicate that the break in the spectral index is only 
created by EBL attenuation. 

Disentangling intrinsic spectral softening from the effects caused by the 
EBL is complicated, and requires a clear understanding of the physical pro- 
cesses that cause the intrinsic softening or hardening (which is very rare) 
between the GeV and TeV energy regions of the source spectrum. 

Multiwavelength observations are required to solve this problem. Such 
observations will provide conclusive tests of non-thermal 7— ray emission 
models, and constitute an important step towards achieving the ultimate 
goal of using the radio-optical- X-ray and low energy 7— ray spectra as a pre- 
dictor for the TeV spectrum. The application of this approach to different 
source classes, i.e., blazars, radio galaxies and starburst galaxies could pro- 
vide additional redundancy to help constrain the EBL. 

The current constraints on the UV/optical to mid-IR regions of the EBL 
have come predominantly from studies of blazars, since the vast majority of 
the extragalactic GeV/TeV 7— ray sources are FSRQs and BL Lacs. AGN 
populat ion studies with Fermi reveal only a small number of non-blazar 
AGNs ( Ackermann et aD . 2011 ). Consequently, we will focus in the follow- 
ing sections on EBL limits that were derived from studies of blazar spectra, 
which have already provided important limits on the near- to mid-IR spectra 
of the EBL. 



3. 3. Phenomenology of Blazar Spectra and Models for TeV r y— ray Production 

The non-thermal emission spectra of blazars generally exhibit two emis- 
sion peaks in z/ F, y , the power emit ted per unit logarithmic photon energy 
[see Figure ((H]), ( Abdo et al. . 2011c! )]. The peak in the radio- UV-X-ray wave- 
band is unequivocally attributed to synchrotron radiation that is produced 
by ultra-relativistic electrons. The second peak, located at X-rays or 7— ray 
energies, is commonly believed to stem from soft photons that were upscat- 
tered by the inverse Compton process to X- and 7— ray energies by the very 
same electrons responsible for the synchrotron emission. This mechanism for 
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Figure 6: The SED of Mrk 421 depicts the synchrotron peak at X-r ay energies, a nd the 
inverse Compton peak at TeV energies. The figure was taken from (jAbdo et all [2011c, 
Figure 11). 



creating the second peak is often referred to as the synchrotron- self- Compton 
(SSC) mechanism. If the second peak includes upscattered photons drawn 
from other ambient radiation fields the mechanism is called the external- 
Compton (EC) mechanism. 

Alternatively, the second emission peak can be produced by a hadronic 
jet containing a significant amount of energy in ultra-high-energy (UHE) 
protons that subsequently interact with soft photons via p + 7 —> vr°, ti ± 1 .... 
interactions. The pions decay giving rise to 7— rays, electrons, muons, and 
neutrinos. The strong magnetic fields required to collimate the hadronic jets 
also lead to the generation of considerable synchrotron radiation by protons 
and charged leptons in the pair cascade. The electrons from the pair cascade 
contribute to the lower energy synchrotron peak, whereas the muons together 
with the secondary photons from neutral pion decay contri bute to the higher 
ener gy 7— ray peak . For a recent review of blazar models see lBoettcherl (120101 ) 
and iDermerl ( 120121 ) . 

BL Lacs are classified loosely by their synchrotron peak position. They 
are referred to as low-frequency peaked BL Lacs (LBLs) if their synchrotron 
peak is at 10 13 < v vea ^ < 10 14 Hz; as intermediate-frequency peaked BL 
Lacs (IBLs) if 10 15 < z/ pca k < 10 16 Hz; and as high-frequency peaked BL Lacs 
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(HBLs) if fpeak > 10 17 Hz (INieppola et al.l . 120061 ). The corresponding 7— ray 
peak follows a similar pattern, with the peak energy progressively shifting 
towards higher energy as the synchrotron peak shifts to higher frequencies. 
LBLs peak at a few GeV (similar to FSQRs), IBLs in the tens of GeV, and 
HBLs generally peak beyond 100 GeV. 

As the 7— ray peak is found at increasingly higher energy, its luminos- 
ity also decreases relative to that of the sy nchrotron pe a k, a trend often 



referred to as the Fossati blazar sequence, (IFossati et al.l . Il998l ). As a re 



suit, the multiwavelength spectra of FSRQs and LBLs exhibit a prominent 
and large 7— ray luminosity, whereas in IBLs and in particularly HBLs, the 
7— ray peak is typically dwarfed by the synchrotron emission. This makes 
both FSRQs and BL Lacs very useful and complementary for EBL studies. 
FSRQs are extremely bright in the GeV regime and despite their low peak 
energy, some are still detectable in the sub-TeV regime with atmospheric 
Cherenkov telescopes. They are thus becoming increasingly useful for EBL 
constraints in the UV/optical/near-IR. HBLs are at the other extreme, their 
7— ray peak can extend well into the multi-TeV regime, while their 7— ray to 
synchrotron luminosity ratio is much smaller. However, during strong flares 
some HBLs have been found to show a 7— ray dominated spectral energy 
distribution. They are therefore useful for constraining the EBL at near- to 
mid-IR wavelengths. 

In general, HBLs can be well described by basic SSC models where the 
ultrarelativistic electrons and their target photons are closely linked via syn- 
chrotron radiation. IBLs are better described by external Compton (EC) 
models, that include a strong ambient photon field external to the blazar jet, 
thereb y providing additional target photons for IC scattering ( Acciari et all 
2009eJ). FSRQs, such as 3C 279, are difficult to fit with either model and 



may require the addition of a hadronic component. 

The detailed fitting of blazar spectra with models requires extensive mul- 
tiwavelength monitoring including optical, X-ray, and 7— rays. Results from 
such multi-wavelength campaigns carry the potential to reduce the spread 
in the Ar GeV (z) relation in figure El While numerous successful multiwav e- 
length campa i gns have been reported in the literature (Abdo et al. , 2011al lcl; 



Aleksic et al.l . l2010al ; IFossati et al.l . 120081 ; lAbramowski et al.l . l201ll ) , the un 



equivocal interpretation of the spectral energy distributions with a clear pre- 
diction for the intrinsic TeV spectrum has not been possible. 

In spite of the fact that blazar spectra are complex, it is possible to set 
limits on the behavior of their GeV/TeV spectra that are based on funda- 
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mental physical limits imposed by energy losses in the particle acceleration 
and radiation processes. The most prominent limitation is the maximum 
hardness of the 7— ray spectrum which provides an important constraint on 
the EBL. The energy spectrum of electrons produced in mod els of diffusive 
shock acceleration in blazar jets (jMalkov &: O'C Druryl . l200ll ). strongly con- 
strains the hardness of the resulting 7— ray spectra produced in SSC and EC 
models, limiting their power law index to values larger than T « 1.5. 

Even spectra that obey this limit are difficult to produce at higher energies 
where the Klein-Nishina effect softens the energy spectra substantially. The 
detection of relatively hard TeV spectra of blazars with redshift ~ 0.1 — 
0.2 therefore came as a surprise, since the absorption corrected spectra are 



already reaching 



(lAharonian et al 



; he T £a 1.5 limit for the min i mal EBL imposed by the IGL 
2006al: iLevenson & Wright) . 120081 : iKrennrich et all 12008 



Ackermann et al. 



20111 : lAbdo et al.l . l2010af ). 



To which extent a problem of hard TeV spectra persists depends on the 
theoretical scenarios i nvoked to explain t h ese sp ectra. A solution to the prob- 
lem was proposed by iKatarzyhski et al.l ( 20061 ). In their model, a high low- 
energy cutoff in the electron distribution could give the app earance of a hard 
7— ray spectrum for a given energy regime. Other ideas by lAharonian et al. 



(120081 ) show that 7 — 7 absorption in the source due to narrow band emis- 
sion from the AGN could lead to unusually hard TeV spectra. T he emis- 



sion produced by proton synchrotron radiation (lAharonian et al.l . 12007a 



Zacharopoulou et al.l . 1201 ll ) combined with internal absorption at lower ener- 



gies has been shown to produce spectra that exceed the hardness achievable 
by diffusive shock acceleration. Other explanations avoid substantial EBL 
absorption by introducing axion-like particles that couple with photons in 
intergalactic magnetic fields thus reducing the 7— ray opacity of the universe 
substantially (jSanchez-Conde et al.l . 120091 ). Mechanisms that would substan- 
tially weaken EBL limits derived from 7— ray observations are discussed in 
Section 8. 



3.4- Unphysical Blazar Spectra 

In spite of the difficulties and the complexity associated with modeling the 
intrinsic 7— ray spectra of blazars, models generally predict a common feature 
of blazar spectra: they follow a power law with curvature and/or exponential 
cutoff and overall can be described by a concave shape (curvature in energy 
flux is downward rather than upward). This means that any absorption 
corrected 7— ray spectrum showing an exponential rise cannot represent a 
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physical source spectrum, and must have its origin in an over-correction for 
EBL absorption.]. 

The relation between the intrinsic, (dN/dE)i nt , and observed, (dN/dE) b s , 
blazar spectrum is given by: 



dEJint 



dEj obs e 



The equation illustrates that an overestimate of the 7 — 7 opacity will lead 
to an exponential rise in the inferred intrinsic spectrum of the blazar. Such 
exponential rise is unphysical. It runs contrary to our basic understanding 
of blazar mode ls, and is absent in the 7— ray spectra of blazars at energies 
below 10 GeV (lAbdo et al.l . l2010af ). This behavior can therefore be used to 
exclude EBL models with optical dep ths that will result in an exponential 
rise i n the corrected blazar spectrum ( iGuy et al.l . l2000l : iDwek fc Krennrichl . 
20051 ). 

Similarly, EBL scenarios can lead to an absorption correction for which 
the reconstructed intrinsic spectra follow a power law with an extremely hard 
spectral slope. Additional constraints can be derived from the spectral slope 
itself, however these are model dependent. EBL models le ading to slopes of 



r < 1. 5 for the intrinsic spectrum, have been rejected by lAharonian et al 



( I2006a[ ) on the basis of the diffuse shock acceleration model. Caveats to this 



approach have been ex 


ensively discussed in the literature ( 


Stecker et al.. 


2007; 


Katarzvhski et al. 


2006: 


Bottcher et al.. 


2008; 


Aharonian et al. 


2008; 


Krennrich et al., 


2008 




Lefa et al. 


> 


2011 


; Zacharopoulou et al.. 


2011 


, and 



while some extreme blazars may exhibit harder spectra with T < 1.5, most 
GeV energy spectra of b lazars obey the T = 1.5 limit, with few exceptions 



( lAckermann et al.l . 1201 if ). 

These general constraints on the hardness of the intrinsic blazar spectra 
have been used to derive c onstraints on the EBL by usin g; few individual ob- 



jects ( IDwek &; Krennrichl. 120051 ; lAharonian et all 
pies of blazars ( Mazin &; Raud . 2007 ; Meyer et al. 



2012 



2006a ) , and by using sa m 



Orr et al.l . 120111 ) 



1 This assumes that no other external effects such as pileups due to pair halos are 
playing a major role in the modification of blazar spectra. The production of pair halos 
would require a magnetic field intensity between f0~ 7 and I0~ 12 G, sufficiently large so 
that electrons are isotropized, and sufficiently small so that IC losses exceed synchrotron 
losses. Instead of a power law with curvature term, a power law with exponential cutoff 
or a broken power law is also consistent with the above statement. 
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4. The Extragalactic Background Light I: Measurements 
and Limits 



4-1- Spectral Measurements and Limits 

The specific intensity of the EBL is usually presented in units of nW m -2 sr" 
The conversion between these units and others that are sometimes used in 
the literature is given by: 

z/4(A)[nWm- 2 sr- 1 ] = -^1 j„(A) [MJy sr" 1 ] (12) 

9.85 x 1(T 3 T ... r , _ 2l 

The EBL intensity can be determined in several ways. The first con- 
sists of direct measurements, a method that poses considerable technical 
and astronomical challenges. Technically, it requires the absolute calibra- 
tion of the instruments, and the understanding and removal of all measure- 
ments uncertainties. Astronomically, it requires the removal of strong fore- 
ground emission from interplanetary dust particles (the zodiacal light, ZL) 
and from stellar and interstellar emission components in the Milky Way. A 
thorough review of the challenges in determining the EBL was presented by 



Hauser fc Dwekl (120011 ) . 

A strict lower limit to the EBL intensity can be obtained by adding the 
light emitted by resolved galaxies. In principle, the integrated galaxy light 
(IGL) can converge to the total intensity of the EBL. A necessary condition 
for convergence is that the spectral index a of the differential galaxy num- 
ber count versus flux S, dN/dS ~ S~ a , becomes smaller than 2 at lower 
fluxes, so that the total integrated intensity, J S 2 (dN/dS) dS, is finite. At 
short wavelengths the intensity of the IGL is limited by the sensitivity of 
the survey. However, even in deep surveys the convergence of the IGL does 
not ensure the measurement of the total EBL intensity, since the low surface 
brig htness regions of galax i es may be missed in stand ard aperture photome- 



try feernstein et all . 120021 : iLevenson fc Wrightl . l2008h . Furthermore, a truly 



diffuse background will always remain undetected in such surveys. 

At longer wavelengths and large beam sizes, unresolved galaxies become 
a source of confusion, limiting the depth of the survey. Below a certain flux 
(the confusion limiting flux) individual sources become indistinguishable from 
the variation in the sky brightness caused by statistical fluctuations in the 
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number of faint resolved or unresolved sources (IDole et al.l . 12004 . and refer- 
ences therein). These limitations can be partially circumvented by stacking 
analysis. Stacking of astronomical images of sources detected at one wave- 
length enhances their signal rel ative to the rando m background fluctuations 
at some other wavelength (e.g. IDole et al.l . 120061 ). The integrated light ob- 



tained by this method is thus closer to the EBL intensity than that obtained 
by integration down to the confusion limit. 

Finally, given a model for dN/dS, one can extrapolate the differential 
source count to very faint fluxes, and evaluate the sensitivity of the integrated 
intensity to the lower flux limit and functional shape of the extrapolation. 

Tables 3-5 list measurements and limits on the EBL intensity derived by 
the different methods described above with the different satellites, balloons, 
and ground observatories. Select measurements were used to define the gray 
area in Figure [7J Absolute measurements and their lcr uncertainties were 
used to define the upper limits on the EBL. The integrated light from resolved 
galaxies and their la uncertainty was used to define the lower limit on the 
EBL. Lower limits derived from stacking analysis were used when available . 



At 140 /im the DIRBE detection with the FIRAS calibration (IHauser et al 
1998), and at longer wavelengths the FIRAS detections by (IFixsen et al. 



19981 ) were used to define the limits on the EBL. The measurements used 
to define the upper and lower limits on the EBL are shown as bold entries 
in Tables 3-5. The figure shows that the EBL is poorly determined in the 
~ 5 — 60 /im wavelength region, whe re the foreground e mission from the 



interplanetary dust cloud is strongest ( IKelsall et al.l . |l998J) 



4-2. Integral Constraints on the EBL Intensity 

The total EBL intensity per unit solid angle is given by: 



l EBL 



C 
47T 



C(z) 



dt 



dz 



dz 
l + z 



(13) 



where C(z) is the luminosity density in a comoving volume element at redshift 
z, and 1 is given by eq. (E])- 

The comoving luminosity density is dominated by the radiative output 
from stars. On a galactic scale, an AGN can dominate the optical to IR 
output of a galaxy, however, on a global scale AGN make only a small con- 
tribution to the total e nergy releases in the un i verse. AGN make up most o f 
the X-ray background (IMushotzky et all I200CH : iDraper &: Ballantvnd. 120 09). 
and a significant fraction of the cosmic radio background (IDwek Barkerl . 
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Figure 7: Limits on the EBL intensity. Lower limits (blue squares) are determined by the 
intensity of the IGL. Upper limits (red circles) are determined by absolute measurements 
of the EBL. The data used in the hgure are listed in Tables 3-5 in bold. The shaded area 
depicts the range of the allowed EBL intensity as determined by UV to sub millimeter 
observations. 



20021 ) . However, they make only a small contribution to the total IR back- 



ground. Accretion onto a central black hole releases about 10% of the rest 
mass energy of the accreted matter, significantly more than the 0.7% re- 
leased in nuclear processes. However, the mass locked up in BH is only 
~ 0.6% of that in stellar objects. Rest frame color-color diagrams generated 
from Spitzer IRAC and MIPS observations covering the 3.6, 4.5, 5.8, 8.0, 
24, and 70 /im bands, of radio-detected submillimeter-selected galaxies with 
spectroscopic redshifts show that AGN constitute a small fr action, between 
13 an d 19%, of the sample dominating its mid-IR spectrum (jHainhne et al 



20091 ). So with these limits in mind, the total comoving luminosity density is 



a direct measure of the cosmic star formation rate (CSFR) at a given redshift. 

Figure M presents the redshift dependence of the CSFR as determined 
from UV-optical emission lines, [O II], [O HI], Lva, Ha, H/3, and from mid- 
IR, submillimeter, and radio ob servations ( jMadau et al.Ul996tlHopkins fc Beacom . 



2006 ; Michalowski et al. . 2010l). The da ta in the figure were taken from the 



compilation of iMichalowski et al.l (120101 ). who used the standard ACDM pa- 
rameters to calculate volume densities and the iKennicuttl (119981 ) law with a 
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Figure 8: A compilation of the cosmic star formation rate as inferred from UV, H a , mid-IR, 
submillimeter, radio, and Lya observations. Lower limits were excluded from the figure. 
The red curves are analytical approximations representing the upper and lower limits to 
the CSFR. 



Salpeter stellar initial mass function (IMF) to convert luminosity densities 
to star formation rates. The two red lines represent a broken power law 
approximation to the upper and lower limits lines of the observations. For 
consistency with the observational determination of the CSFR, the CSFR 
was converted to a bolometric intensity using the Salpeter IMF and a star- 



burst age of 100 Myr: L M = 7.5 x 10 9 ^ M Q yr" 1 flDwek et all . 120111 ). The 



integrated intensity of the CSFR is then bounded by: 

Iebl = 21-66 nW nT 2 sr" 1 (14) 

Table 6 presents the current limits on the EBL intensity, separated into 
its stellar and dust components, and compares them to this predicted by the 
various EBL models presented in §5. 

4-3. Constraints on the EBL from Fluctuation Measurements 

Most of the EBL is generated by discrete galactic or primordial stellar 
sources. Fluctuations in their number and their clustering properties will give 
rise to spatial fluctuations in the EBL intensity. Studies of the op tical region 



of the EBL via the fluctuations method were first conducted by (IShectman 
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19731. Il974h. and of th e near-IR region using the COBE/DIRBE data by 



Kashlinsky et al.l (119961 ) . The fluctuations do not provide a direct measure- 
ment of the EBL intensity. The derivation of the EBL intensity from these 
measurements will require detailed knowledge of the galaxy source counts, 
their luminosity function, and their clustering properties as a function of 
redshift. However, spatial fluctuations in the EBL provide a different means 
of setting limits on its intensity. Fluctuation measurements of the EBL do 
not require absolute measurements of its intensity, since the removal of fore- 
ground emission components is done on the basis of their distinct spatial 
properties rather than their absolute intensities. At IR wavelengths, these 
fluctuations can have a spatial signal that is distinctly different from that 
generated by the interplanetary dust cloud or by interstellar dust. After the 
removal of all known resolved sources, any residual fluctuations will measure 
the EBL contribution from sources that may represent a yet unknown or 
unresolved population of stars or gal axies. Recent fluc t uatio n measurements 



in the Spitzer 3.6 and 4.5 /im bands (iKashlinsky et al.l . 120121 ) show that they 



are in excess of those that can be attributed to known galaxy populations 
(jHelgason et al.l . 120121 ) . suggesting an origin in a very faint, yet unknown, 



population of highly clustered sources. More details and references on the 
use of fluctuations ana lysis to constrain the E BL intensity can be found in 



Hauser fc Dwekl (l200ll ) and iKashlinskvl (j2005j). 



EBL fluctuations at far-IR w avelengths were detec ted at 170 /im ( jLagache fc Puget 



2000 ', ISO data) 



at 160 /zm (jLagache et al.l . 120071 . Spitzer data) and by 



Shang et al.l ( 120121 ). using th e Planck dat a . Flu ctuation measurements and 



EBL colors have been used by lPenin et al.l ( 120121 ) to derive limits on the EBL 
at 100 and 160 /im (see Tables 4 and 5). 



5. The Extragalactic background Light II: Models 

The EBL intensity only provides an integral constraint on all the radiative 
energy releases over cosmic time. In a dust-free universe, it represents the 
intrinsic stellar or AGN spectra. The EBL intensity and spectral shape de- 
pends then on the star formation history, the stellar initial mass function, the 
evolution of metallicity, the energy released by AGN, and the relative impor- 
tance of energy releases by nuclear and gravitational processes. In a dusty 
universe, the comoving luminosity density, C(z), and total EBL intensity, 
I ebl remain unchanged, however the energy is redistributed by absorption 
and reemission processes over a large spectral range. The resulting spectrum 
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depends on many factors, ranging from the size distribution, composition, 
and optical properties of the dust grains, the evolution of their abundance 
and properties over time, and on the morphology of the galaxy which deter- 
mines the spatial distribution of the dust relative to the radiative sources. 

Several distinct approaches have been used to model the intensity and 
spectral distribution of the EBL at z — 0. They all represent different ap- 
proaches for calculating the evolution of C v (X,z) with redshift [see eq. (2)]. 
Backward evolution models start from the local determination of C u (X,z = 
0), evolving it with redshift using observed galaxy number counts at different 
wavelengths. Forward evolution models use the CSFR to determine £„(X, z) 
as a function of redshift, and population synthesis and radiative transfer 
models to determine the distribution of the energy over wavelengths. Cosmic 
chemical evolution models are similar to previous models, except that their 
system is the universe as a whole. Semi-analytic models calculate C V (X, z) 
by including the appropriate physical processes in a more general model 
for the formation and evolution of structure in the universe. In the follow- 
ing we discuss these different models in somewhat more detail. For a more 
extensive discus s ion a nd description we refer the reader to the review by 



Hauser & Dwekl (120011 ). 



5.1. Backward Evolution (BE) Models 

BE models start with the construction of a library of the galactic spec- 
tral energy distributions (SEDs) representing those of galaxies in the lo- 
cal universe and evolve them back in time in order to fit observed number 
counts. The EBL intensity and spectral shape is used as an integral con- 
straints on their evolution. The galaxies in such library should represent the 
range of observed galactic morphologies (spiral elliptical, irregular) and ac- 
tivities (AGN, normal, starburst, mergers) in the local universe. The SEDs 
of the galaxies comprising such library should also satisfy observed statistical 
properties of the ensemble of local galaxies, such as: the trend of increas- 
ing S (60 urn) / S (100 jj<m) /) and decreasing S (12 p,m)/S (25 ur n)/) flux ratios 



with increasing IR luminosities (( jSoifer fc Neugebauerl . Il99ll ); and the num 



ber density of galaxies in the L + dL luminosity interval, represented by the 
luminosity function (LF), <&(L)dL. 

Many functional forms have been adopted to characterize the local LF. 
The most commonly one used at optical and near-IR wavelengths is the 
Schechter (1976) LF. At wavelengths above ~ 10 /xm the galaxies' SED is 
dominated by thermal emission from dust, and their LF cannot be adequately 
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represented by the Schechter LF. Several distinct functional forms have been 
used to characterize the LF in the different IR wavebands. They are usually 
characterized by three parameters: a normalization parameter, a char- 
acteristic luminosity, L*, that determined the transition point between the 
low and high luminosity behavior of the L; and a power law index, a that 
determines the behavior of the LF at low luminosities. A list of references to 
the functional forms derived from the IRA S s urvey can be found in Hauser & 
Dwek (2001). Of those, the parameters of the lSaunders et al.l ( 119901 ) LF have 
been recently updated t o fit the differential 24 p ,m number counts obtained 
by deep Spitzer surveys (IRodighiero et al.l . 120101 ) . 

If neither the galaxies' SED nor their comoving number density evolved 
with time, the spectral luminosity density, £ v (\, z) would be independent 
of redshift. However, the recent deep surveys with the Spitzer and Her- 
schel satellites show strong evolution in number counts, compared to predic- 



tions made with no evolution models ( ILagache et al.l . l2005t IRodighiero et al 



2010l ). Evolution in the LF or, equivalently, the spectral luminosity density 



can be inferred directly from observations if the r edshift of the sou r ces is 
known, and t heir n umber counts are complete (e.g. IRodighiero et al.l . 12010 ; 



Dunne et al.l . 120001 ) . Alternatively, evolution in the LF can be introduced 
by adding a redshift dependence, usually of the form (1 + z) 1 , to the basic 
parameters, L+, and a, that characterize the LF. The value of 7 is then 
derived by fitting model prediction to the observed galaxy number counts in 
a given waveband. Evolution in the LF can also be modeled by evolving the 
relative number of the different type galaxies: quies cent star forming galax- 



ies, s t arbursts, AGN s , and ellipticals with re d shift ( IRowan- Robinson! 12001 



20091 : ILagache et al.l . 120031 : iDommguez et al.l . l201ll ). or by simply evolving 



the relative number of ultraluminous infrared galaxies (ULIRGs), character- 
ized by IR luminositie s in excess of ~ 10 12 L^, relative to the rest of the 



galaxy population [e.g. IChary fc Elbazl ( 120011 )]. Having determined the evo- 
lution of the spectral luminosity density with redshift, the EBL is obtained 
by a simple integration of CJ\ z ) over redshift. 

Hekason fc Kashlinskvl (120121 ) adopted a more empirical approach than 



standard BE models, using number-count derived luminosity functions at 
UV (0.1 fim) to mid-IR (25 jum) bands, to construct the luminosity density, 
£„(\, z), as a function of redshift. 

BE models are relatively simple, and their predictions can be easily com- 
pared to observed galaxy number counts, their magnitude-color and magni- 
tude number density relations, and their redshift distribution. They are only 
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loosely constrained by physical processes. The predicted dependence of the 
comoving bolometric luminosity density with redshift should be consistent 
with that inferred from limits and observations of the cosmic star formation 
rate. 



5.2. Forward Evolution (FE) Models 

FE models use the redshift dependence of the cosmic star formation 
rate (CSFR), inferred from a variety of way ebands and line observation 
( jMadau et all 1 19961 ; iHopkins fc Beacoml . 120061 ) . as a starting point in their 
calculations. The CSFR is determined from the UV-optical line and con- 
tinuum, and IR and radio emission in the various wavebands using statis- 
tically dj^tejmh^d galaxies in the local uni- 
verse flKennicuttl . Il998l : lHaarsma etUl l2000h . Determination of the CSFR 



is complicated by extinction effects at UV and optical wavelengths, and by 
the implicit assumption that the IR luminosity is powered by stars and repre- 
sentative of the total bolometric luminosity of the galaxies. Even if the total 
bolometric luminosity of a given galaxy is determined, its conversion to a star 
formation rate requires knowledge of the stellar IMF, a poorly determined 
quantity at high redshifts, and the duration of the starburst activity at each 
redshift. Once the CSFR is determined, FE models use popu lation synthe- 
sis models such as PEGASElFioc fc Rocca-Volmerangel (119971). Starburst99 



Leitherer et al.l (119991 ) , iBruzual fc Charlotl (120031 ) or iMaraston fc Stromback 



( 1201 it ) to calculate the stellar bolometric and spectral luminosity density as 
a function of redshift. 

The most difficult part of this approach is determining the fraction of 
starlight that is absorbed by dust, and the spectrum of the reradiated IR 
emission. The SED of a galaxy can be det ermined w ith radiative transfe r 
models, suc h as GRASIL dSilva etldl Il998h. DIRTY dGordon et aD . l200lh . 



or DUSTY ( Nenkova et al. . 2000l ). Alternatively, one can use a parametric 



approach in which the fraction of UV-optical light absorbed by the dust, and 
the reradiated infrared spectrum are statistically determined from observa- 
tions. EBL spect r a derived from FE evo lutio n models we r e pre sented by 



Dwek et al.l (Il998f ). iRazzaque et al.l ( 120091 ). and lFinke et al.l (l2010f ). 



Population synthesis models, combined with simple radiative transfer cal- 
culations, are useful for determining the UV to radio SED of individual galax- 
ies. However, the cosmological application of such models assumes that star 
formation is a monolithic process, in which in all galaxies star formation com- 
menced at the same redshift and evolved quiescently until the current epoch. 
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The models do not allow for galaxy interactions, stochastic star formation 
histories associated with merger events, or any morphological evolution of 
galaxies. Any discrepancies between model predictions and galaxy number 
counts must be introduced in an ad hoc fashion by evolving the stellar IMF, 
or by introducing a new population of galaxies at the appropriate redshifts. 

5.3. Cosmic Chemical Evolution (CCE) Models 

CCE models treat the universe as a closed system in which all galax- 
ies within a large comoving volume element are represented by their basic 
ingredients: stars, interstellar gas, metallicity, and radiation. Chemical evo- 
lution equations, a nalogous to those us e d to f ollow th e chemical evolutio n 
of the Galaxy (e.g. lAudouze k Tinslevl . Il976l : iTinslevl . Il98ll : IPagell . bOQlf ). 
are used to follow the evolution of the average stellar, gaseous, and radia- 
tive contents in each com oving volume in a self consistent manner. CCE 
models w e re pio neered by iPei fc Fall! (119951 ). and most recently updated by 
Pei et al.l (119991 ). Inputs parameters for their model are the mean rest frame 
UV luminosity density as a function of redshift, and the mass of the ISM 
gas as determined from H I column densities derived from studies of quasar 
absorption lines through damped Lya systems. The decrease in the ISM gas 
with redshift and the UV luminosity density were used to derive a solution 
for the evolution of the CSFR with redshift which is consistent with that de- 
termined from the extinction-corrected Ha, and with SCUBA 850 and ISO 
15 /im surveys. Similar to FE models, population synthesis models were 
then used to calculate the stellar SED at each redshift, and an LMC extinc- 
tion law was adopted to calculate the fraction of starlight absorbed by the 
dust. A power-law distribution in dust temperature was used to calculate 
the spectrum of the reradiated IR emission. The model reproduced various 
observational constraints, including the comoving rest-frame 0.44, 1.0, and 
2.2 fim spectral luminosity densities in the ~ — 2 redshift interval, the 12, 
25, 60, and 100 /im local luminosity densities; and the mean abundance of 
metals in damped Lya systems in the ~ 0.4 — 3.5 redshift interval. 



5.4- Semi-analytical (SA) Models 

SA models follow the formation and evolution of galaxies in a cold dark 
matter Lambda dominated (ACDM) universe using the cosmological parame- 
ters derived fr om the 5-year Wilkins on Microwave Anisotropy Probe ( WMAP5) 



observations (IHinshaw et al.l . 120091 ) as the initial conditions. SA models then 



follow the growth and merger of dark matter halos, and the emergence of 
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galaxies which form as baryonic matter falls into the potential wells of these 
halos. The fate of the infalling gas is determined by many different pro- 
cesses: the formation of stars in a multiphase interstellar medium, AGN and 
supernovae feedback processes that quench their formation, the evolution of 
the stellar radiation field, the heating and cooling of the interstellar medium 
and its chemical enrichment, the exchange of material with the intergalactic 
medium through infall and galactic winds, and the growth of the central black 
hole. A descripti on of recent deve l opmen ts and references to previous work 



can be found in ISomerville et al.l ( 120111 ). Model prediction are compared 
to a basic set of observational constraints such as the observed characteris- 
tics of galaxies: their morphology, colors, and spectral energy distribution, 
and morphology; and their integrated cosmological properties: their number 
counts and luminosity function in different wavebands and redshifts, their 
mass function, the cosmic star formation rate, and the EBL generated by 
them. As in all EBL models, determination of the galaxies' SED is com- 
plicated by the detailed microscopic and large scale parameters needed to 
calculate the amount of starlight that is absorbed by dust, and the spectrum 
of the reradiated emission. Recent SA models have combined the models for 
galax y formation with radi a tive transfer models to determine the galaxies ' 



SEP (iFontanot et all 120091; IFontanot fc Somervilld . 1201 It ISomerville et al. 



201ll : lYouneer fc Hopkinsl . l201lh . 



SA models are inherently complex, incorporating a large number of phys- 
ical processes, some poorly known, to derive galaxy properties. However, 
they are the most physically motivated models, and quite successful in re- 
producing a large number of observational constraints. 



5. 5. Comparison of Model Predictions with Observations 

A detailed comparison of all model ty p es wi th EBL limits and obser- 
vations was presented by lHauser fc Dwekl (120011 ). Here we will represent 
mostly the mod e ls tha t have been developed s ince then: BE models by 



Stecker fc Scully! (120061). iFranceschini et al.l ( 120081 ) , and lDommguez et al.l (120111) : 



The F E model of iFinke et al.l ( 120101 ): and the SA model of iGilmore et al. 



( 1201 it ). Figure IH] compares the various models to the current limits and ob- 
servations of the EBL. In general, all models, except for the BE models of 
Stecker et al. provide adequate fits to the EBL. 
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Figure 9: Models of the EBL are compared to observational limits on the EBL. 



6. EBL Constraints from 7— ray Observations of Blazars 



The attenuation of 7— rays by the EBL can in principle be used to deter- 
mine the EBL intensity at wavelengths corresponding to the 7— ray observa- 
tions. Neglecting the possible scattering or production of second generation 
7— ray photons along the line of sight to the blazar, the intrinsic 7— ray flux 
from the blazar, F Int (E 7 ), can be related to the observed one, F obs (E^) by: 



F o6s (E 7 ) = F Int (Ery) exp[-r 77 (E 7 )] 



(15) 



where the optical depth, r 77 , is given by eq. (jUJ). Determination of the EBL 
assumes that all the attenuation is caused by interaction with the EBL, in- 
stead of photons in or around the vicinity of the blazar. Furthermore, it 
requires knowledge of the intrinsic blazar spectrum. Assuming that all the 
attenuation is attributed to the EBL, several upper limits have been derived 
on the EBL intensity by making various assumptions on the intrinsic blazar 
spectrum. The 7— ray derived EBL limits are compared to those derived 
from UV to sub millimeter observations in Figure [101 and described below. 
Fixed power law: Early observations of the blazars Mrk 421 and 3C 279 
suggested that their ~ GeV-TeV spectrum could be ap proximated by a sin- 



gle power law ( IPunch et all Il992t IStecker et all . Il992l respectively). If so 



then any deviations of the observations from the extrapolated power law to 
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higher energies should be attributed to EBL attenuation. IStecker fc de Jager 
(119931 ) derived an upper limits at 1-5 /im of 10 nW m -2 sr -1 , assuming 
that a straight power law extrapolation of the spectrum of Mrk 421 from 
GeV energies obtained from the EGRET fjLin et all Il992l ) with an index of 
T = 1.96 ± 0.14 holds up to the TeV regime, where t he index was measu red 
by the Whipple collaboration to be T = 2.25 ± 0.19. biller et all (Il995h in- 
cluded the statistical uncertainties of the GeV spectrum, and demonstrated 
vastly different extrapolations with significantly higher upper limits, thereby 
yielding conservative upper limit to the EBL in the mid-IR at 10/im. 
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Figure 10: Limits on the EBL as determined from 7— ray observations of blazars. Details 
in text. 



Strong flaring activity of Mrk 501 provided a well measured TeV spec 



trum from 0.2 - 24 TeV (lAharonian et all Il999l; ISamuelson et all 11998 ; 



Djannati-Atai et all 119991 ). IStanev &: Franceschini (1998) calculated limits 



by fitting the energy spectrum of Mrk 501, with a range of possible absorp- 
tion scenarios assuming intrinsic power law spectra and varying levels of EBL 
intensity by scaling the lower limits f rom galaxy c ounts . This provided strong 
EBL limits in the near- and mid-IR. iFunk et al.l (|1998l) followed a simi l ar ap- 
proach, except that they used an EBL model by lMacMinn fc Primackl (119961 ) 
as the basis for scalin g the EBL intensity, yielding similar results in the mid- 
IR. iMannheiml (119981 ) argued that the observed spectra would deviate from a 
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power law, if the primary 7— ray spectrum were substantially attenuated by 
the EBL. An upper limit in the mid-IR, based on this hypothesis and on an 
energy spectrum of Mrk 501 from the H EGRA collaboration, is also shown 
in Figure [TOj Later on, IVassilievI ( 120001 ) demonstrated that the absence of 
deviations from a power law does not preclude the presence of substantial 
absorption in the observed spectra. 

The obvious drawback of the method is that the assumption of a single 
power law for the intrinsic blazar spectrum does not hold true over a wide 
range in energy. Blazar spectra generally exhibit a concave 7— ray peak over 
a sufficiently large energy range. Most blazar spectra measured by Fermi or 
by Cherenkov telescopes can be represented by a power law over the energy 
range covered by the instrument. 

Synchrotron Self-Compton (SSC) Spectrum: The SSC model is a pop- 
ular model explaining the existence of the two peaks in the blazar spectrum: 
the synchrotron peak at radio-UV-X-ray energies and the inverse Compton 
(IC) peak at 7— ray energies. The spectrum of the IC peak can be modeled 
using parameters that produce the synchrotron peak and th e unab sorbed 
part (E < 10 GeV) of the IC spectrum (see review bv lDermerl . 120 12h . 



Such models for the intrinsic blazar spectrum have been used by lGuy et al. 



( 120001 ) to determine the intensity of the EBL in the 1-5 /im and 20-80 /im 
wavelength region. They applied a multiwavelength fit to the X-ray and TeV 
data of Mrk 501 in the framework of a standard homogeneous SSC model 
to derive the level of absorption present in the TeV spectrum. As a result, 
they obtained an absolute upper limit on the EBL of 60 nW m -2 sr _1 and a 
most likely value of 20 nW m -2 sr" 1 at 1 fim. They also pointed out that the 
lack of an absorption signature in the spectrum of Mrk 501, as suggested by 
the HEGRA telescopes, does not necessarily imply a lack of EBL absorption. 
They emphasized that in the transition region from the near-IR to the mid-IR 
EBL, the opacity could be nearly constant. This is a consequence of the large 
width of 7 — 7 cross section (see Figure 2). So when <r 77 is convolved with 
the number density of background photons, any strong wavelength variations 
in the EBL are smoothed out. As a result, the observed TeV spectrum at 
1-10 TeV would corresponds to the intrinsic blazar spectrum since the ob- 
served spectrum is now described by (dN/dE) int x e T with r a slowly varying 
function of energy. 

The drawback of using 7— ray emission models to constrain the EBL is the 
uncertainty in the many parameters that determine the IC spectrum. Fur- 
thermore, while HBLs generally can be well fit by SSC models, IBLs require 
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the inclusion of additional ambient radiation fields that make a contribution 
to the 7— ray IC component. 

Additional complications arise from the fact that basic one-zone SSC 
models are not applicable for sources exhibiting "orphan flares", where only 
the TeV flux is enhanced while the synchrotron emission remains unchanged 
( IKrawczynski et al.l . 120041 ). Finally, the biggest challenge for the SSC/multi- 
wavelength approach to constraining the EBL is to get simultaneous mea- 
surements for large sets of blazars. 

The T > 1.5 limit on the hardness of the blazar spectrum: A more re- 
laxed assumption on the intrinsic blazar spectrum is that it cannot produce 
too many hard photons, so that the 7— ray spe ctrum, expressed as dNldE y 
E~ r cannot be flatter than one with T = 1.5 dMalkov fe O'C Drurvl . bpOlh . 
In the spirit of this limit to the spectral index, iRenault et al.l ( 1200 ll ) explored 
a range of EBL scenarios based on measurements with the minimal assump- 
tions that the intrinsic power of Mrk 501 is concave, effectively requiring a 
decreasing energy flux distribution above 4 TeV (T > 2.0). They derived an 
upper limit of 5 nW/m 2 /sr at 10 fim. 

The strict assumption of T > 1.5 , was used by lAharonian et al.l (j2006al ) to 
derive upper limits on the 1-5 /im on the EBL which are close to the lower 
limits determined by the IGL, suggesting that the EBL h as been largely 
resolved at these wavelengths. A comprehensive study by iMazin fc Raue 
( 120071 ) is based on eleven blazars over a redshift range from 0.03 - 0.18, and 
explores a large number (8 million) of hypothetical EBL scenarios to set 
upper limits on the EBL, again with the requirement that the source spectra 
cannot be harder than T = 1.5 or T = 2/3. The lower value arises from the 
extreme scenario of a mono-energetic energy distribution of ultra-relativistic 
electrons in which the resulting IC 7— ray spectrum could be as hard as 
7 = 2/3, leading to two conditional u pper limits. The first con dition yielded 
limits that are slightly above that of lAharonian et al.l ( 2006al ). The second, 
more relaxed condition, yielded limits that were higher by about 30%. 

The theoretical validity of a strict hardness limit of T > 1.5 has been dis- 



cussed by a number of authors, with no unanimous verdict (iKatarzvriski et al. 



2006 



2011 



Stecker et all.l2007l:[Bottcher et all . l2008l : lAIiaronian et"ai" . 20081: Lefa et al 



Zacharopoulou et al.l . l201ll ). Observational evidence, e.g. iLevenson fc Wright 



(120081 ). have provided lower E BL limits from gal a xy co unts that are higher 
than previous ones derived bv lMadau fc Pozzettil (l2000h . If these new limits 
are correct, they imply 7— ray spectra that are slightly harder than T = 1.5 
dKrennrich et all 120081 ) . 
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Unphysical exponential rise of the blazar spectrum: Less model de- 
pendent, and therefore more robust limits on the hardness of the intrinsic 
blazar spectra arise from the notion that an exponential increase of their 
luminosity with energy is unphysical. All current blazar models produce a 
concave spectrum, rendering intrinsic blazars spectra with an exponential rise 
in energy flux theoreticall y unfeasible. The paradigm of concave intrinsic en- 
ergy spectra was used by Dwek fc Krennrich J 20051) to r eject many different 
realization of the EBL. Furthermore, iDwek et al.l ( ]2005bl ) r uled out the extra- 



galact ic origin of the near-IR sky brightness observed by iMatsumoto et al. 
(120051 ). since it would lead to an exponential rise in the spectrum of the 



blazar PKS 2155-304, which is ruled out by observations ( lAharonian et al. 



2005bj). An EBL spectrum close to the IGL limits yielded a blazar spectrum 



consistent with the SSC model, suggesting th at the EBL was mostly resolved 



at near-IR wavelengths ( Dwek et al. . 2005bl JaT). 



Spectral break analysis due to EBL spectrum: lOrr et al.l (1201 if ) de 



veloped a novel approach to set limits on the EBL intensity by studying 
the effects of the spectral shape of the EBL on the sub- to multi-TeV spec- 
tra of blazars. The 7— ray opacity around 1 TeV is sensitive to the EBL 
intensity at ~ 1 /im (see Fig. 3). Its subsequent energy dependence hinges 
on the rate at which the stellar UV, optical, and near-IR emission decreases 
towards mid-IR wavelengths. The near- to mid-IR intensity ratio of the EBL 
determines the relative ~ 1 to 10 TeV opacity. This is illustrated in Figure H] 
which shows the energy dependence of r 77 . A large near- to mid-IR ratio in 
the EBL intensity would cause the 7— ray opacity to be relatively flat in the 
~ 1 — 10 TeV region, resulting in a hard 7— ray spectrum. Conversely, a low 
near- to mid-IR ratio would result in an increase in the ~ 10 TeV opacity 
relative to that at ~ 1 TeV, resulting in a softer blazar spectrum. 

The spectral shape of the EBL in the ~ 1 — 15 /im region can be related 
to the break AT TeV = T(< 1 TeV) — T(> 1 TeV) in the blazar spectra. A 
study of 12 blazars spanning a redshift range from 0.03 to 0.186, showed a 
trend of increasin g | ATxeV I with redshift with a statistical significance of 3.6<r 



( lOrr et al.l . 1201 ll ). This strongly suggests that the trend is caused by EBL 



absorption, providing strong constraint on the near-IR to mid-IR EBL inten- 



sity ra tio. Combined with 7— ray-derived upper limits on the EBL. lQrr et al 



( 1201 ll ) derived correlated constraints between the near-IR and mid-IR inten- 
sity ratio of the EBL and the EBL intensity at mid-IR wavelengths. 
7— ray inverse Compton emission from radio l obes: A new method fo r 



determining the local EBL intensity was described by ( iGeorganopoulos et al. 
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20081 ). The method relies on the detection of 7— rays produced by IC scat- 



tering of CMB and EBL photons off the relativistic electrons of the lobes of 
radio galaxies. Since the lobes have to be clearly resolved with the 7— ray 
telescope, the method is presently limited to the nearby radio galaxy For- 
nax A. With the normalization and maximum electron energy in the lobes 
determined by their synchrotron spectra (measured by WMAP), the 7— ray 
spectrum of the lobes is determined by the intensity of the radiation field 
in the radio lobes, which is dominated by the CMB and EBL. The result- 
ing 7— ray spectrum comprises distinct contributions by the CMB and EBL 
photons: The CMB contribution peaks at energies ~ 40 MeV, while that of 
the EBL appears as an excess above steeply dropping CMB contribution at 
higher energies. This excess emission is relatively flat, and extends to ener- 
gies of ~ 50 GeV. Its magnitude provides a direct measurement of the EBL 
intensity at mid- and far-IR wavelengths. 

Combined 7— ray limits: Figure [11] shows select limits on the EBL de- 
rived from TeV observations. All limits were adjusted to a common Hubble 
constant of 70 km s _1 Mpc -1 . Strict convergence of the different limits is 
not expected, since each limit was derived from different assumptions on the 
intrinsic blazar spectra and the EBL spectrum. The bold horizontal lines 
and the lOrr et al.l (120111 ) result represent the most recent constraints on the 
EBL. It is important to emphasize that different methods were used, i.e., the 
r > 1.5 limit on the hardness of the blazar spectrum and the spectral break 
analysis due to EBL spectrum, yet they reached similar conclusions. 

In general, much improved upper limits are now available from the ob- 
servations with the new generation of atmospheric Cherenkov telescopes 
(H.E.S.S., MAGIC, VERITAS). Both the near-IR and the mid-IR intensity 
levels are now constrained to much lower values than was possible with early 
results from the first few TeV blazars, and make only minimal assumptions on 
the intrinsic blazar spec tra. Extreme EBL scenar ios, with 70 nW m -2 sr _1 at 
1.5 jttm, as suggested by lMatsumoto et al.l (120051 ) are clearly ruled out. Min- 
imal assumptions on the blazar spectra such as the absence of exponential 
rises or applying the r > 1.5 limit, reject such high EBL intensities in the 



ne ar-IR (IDwek et al.l. l2005bl : lAharonian et al.l . l2006al ) 



Mazin fc Rauel (120071 ) provide upper limits that are comparable to the 



former two in the near-IR, while in the mid-IR their upper limits are much 
higher. For this work the spectrum of 1ES 0229+200, a blazar with a redshift 
of z=0.129 and a spectrum up to 10 TeV was not available, while it is the basis 



for strong mid-IR limits in lAharonian et al.l (l2007bl ) and lOrr et al.l (120111 ). 
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Figure 11: The EBL constraints as derived from 7— ray observations to the EBL are shown 
for 1 /im (left) and for 10 /im (right). The horizontal line represents the intensity of the in- 
tegrated galaxy light (IGL), which provides a strict lower limit on the EBL intensity. The 
most recent limits are represented by bold lines. The shaded area re presents the la un- 
certai nty in the IGL intensity. References to t he 7— ray lim i ts are : (0) IStecker fc de Jager 
(Il993l): (Tl IStanev fe Franceschinil (Il998l): (2) iBiller et al.1 (Il998l): (3) iGuv et all (l2000h : 
(A) iDwek et all (l2005bh : (5) ISchroedter et al l (l2005h : (6) lAharonian et al.1 (l2006ah: (7) 



Mazin fc Rauef(l2007h: f S'llFinke fc Razzaaud (l20n9|W9'l lOrr et alJ ^Ollh: (10) 
(Il998h: milMannheiml(ll998t ): Q2^ lDwek fc Slavinl (|l994 : QS'l lR.enault et al 



Aharonian et al.l (12007. 



Funk et all 



20011); (14) 



This emphasizes the importance of extending the TeV energy spectra of 
distant (z> 0.1) blazars into the multi-TeV regime, where 7— rays reach their 
maximum cross-section with photons in the mid-IR. Furthermore, energy 
spectra between 100 GeV and 10 TeV are sensitive to the spectral shape of 
the EBL, thereby linking the upper limits in the near- and mid-IR; a given 
mid-IR intensity level combined with the EBL spectral shape limits the range 
of near-IR EBL intensities, i.e., by excluding EBL intensities that are either 
too high or too low in the near-IR based on the spectral shape constraint. 

There are now in the order of 30 extragalactic objects with redshifts up 
to z w 0.5 available for constraining the EBL. Limits derived from the recent 
generation of experiments ( H.E.S.S., MAGIC, VERIT AS) suggest a l o w mid - 
IR, such as the results from lAharonian et al. I d2007bh and lOrr et all fepllh . 

In summary, it should also be noted that all of the most recent TeV 
constraints are well within the boundaries set by direct measurements and 
their uncertainties (see shaded area in Figure [TO]) . However, given the debate 
about some of the assumptions about the blazar spectra, namely, the T > 1.5 
limit on the hardness of the blazar spectrum, claims that the EBL has been 
resolved are premature. 



35 



7. The 7— ray Opacity of the Universe 



With blazars being detected at increasingly large redshifts, it becomes 
possible to use them to discriminate between different EBL models. The 
7— ray opacity to a blazar at redshift z, is given by: 



T 77 (i£ 7 , z) — — lo£ 



7 Zjobs 

F(E. 



7 jint 



(16) 



The opacity can be determined from models of the EBL if its evolution with 
redshift is known and, independently, from 7— ray observations if the intrin- 
sic blazar spectrum is known. Concordance between these two independent 
determinations of r 77 can serve as a test for the validity of the underlying 
assumptions in each method. 

Figures [12] and [13] depict the evolution of the comoving intensity of the 
EBL, the corresponding evolution of the proper number density of back- 
ground photons, the optical depth to blazars at various redshift, and the 
correspo nding attenuation f actor . Results are plotted for the B E evolution 
model of lFranceschini et al.l (120081 ) and the BE evolution model of lDommguez et al 
fl201lh . 

Determining the 7— ray opacity from observations requires knowledge of 
the intrinsic blazar spectrum. Differences between the observed and expected 
flux at a given energy E 1 would then be simply attributed to EBL attenu- 
ation. Figures [T2] and [13] show that the sharp drop of the EBL intensity at 
UV and shorter wavelengths renders the universe almost transparent to GeV 
photons. Consequently, the observed ~ 1 — 50 GeV spectrum is very likely 
the intrinsic blazar spectrum. So instead of assuming a theoretical limit on 
the spectral index, one can use the GeV - 10s of GeV energy spectral slope 
from Fermi data as a proxy for the intrinsic spectra at TeV energies. 

Assuming that this power law can be extrapolated from GeV to TeV en- 
ergies, one can derive t he TeV optical depth t o the observed blazar. This 
approach was u sed by iGeorganopoulos et al.l (120101 ) and in method 1 in 
Orr et al.l (1201 ll ) to set firm upper limits on EBL models using the GeV 
to TeV spectra of PKS 2155-304 (z = 0.116) and 1ES 1218+304 (z = 
0.182). Assuming that the GeV spectrum is unattenuated by the EBL, 
Mankuzhiyil et al.l ( I2010l ) used optical, X-ray and GeV data to model the 
TeV flux of PKS 2155-304 using a one-zone SSC model. Comparison of the 
model results with observations, they derived the TeV opacity to this blazar, 
and found it to be consistent with most EBL models. 
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Figure 12: Basic EBL model results bv lFranceschini et al. ( 20081) : Top left: The comoving 
EBL and CMB intensities versus wavelength for different redshifts; Top right: The proper 
number density of EBL and CMB photons versus energy for the same grid of redshifts 
as the previous panel; Bottom left: The 7— ray opacity versus energy, E 1 for different 
redshifts; Bottom right: The amount of attenuation versus energy for the same grid of 
redshifts as the previous panel. The figure illustrates the change in the slope of r 77 at 
energies corresponding to the wavelength at which the slope of the EBL spectrum changes. 
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Figure HH compares the dependence of the optical depth derived from 
EBL models (hatched curves) to that derived for select blazars: Mrk 501, 
1ES 1218+304, and 3C 66A. Each ha tched band spans the range of optical 
depth s predicted by the EBL mo dels of iFranceschini et al.l (120081 ) , iFinke et al 
fl2010h . lDommguez et all (boilh . and lGilmore et all feoilh . The colored dots 
represent the optical depths derived from the 7— ray observations of the three 
blazars. The intrinsic blazar spectrum was assumed to be a power law de- 
termined by the observed flux at 1 GeV and the spectral index, Tcev- The 
observed flux in the TeV range was assumed to be a power law with a spec- 
tral index I+ev (see Table 2). The 7— ray opacity in the TeV range was then 
derived from eq. (9). The band of opacities for each blazar was obtained by 
performing 100 Monte Carlo simulations of the intrinsic and observed spec- 
tra using the uncertainties in the spectral indices and 7— ray energies into 
account. 

The figure shows that the 7— ray derived optical depths of Mrk 501 and 
1ES 1218-304 are in general agreement with model prediction. The discrep- 
ancy between the EBL and the 7— ray derived optical depth for 3C 66A is 
typical of most blazars liste d in Table 2. We no te that the redshift to 3C 66A 
is still somewhat uncertain ( Abdo et al. . 2011bl ldl). The convergence between 
observational limits on the EBL and models suggests that the origin of the 
discrepancy can be mostly attributed to our still incomplete knowledge of 
the intrinsic spectra of blazars. 

The EBL not only affects the 7— ray spectra of individual 7— ray sources, 
but also the spectrum of the extragalactic 7— ray background (EGRB) which 
consists of the cumulative contribution of resolved and unresolved sources 
and a possible truly diffuse emission component. 

Recently, the Fermi Large Area Telescope (LAT) provided a new mea- 
surement of t he diffuse 7— ray ba ckground (DGB) at energies between 0.2 



and 100 GeV (lAbdo et al.l . l2010el ). obtained by the subtraction of resolved 



Fermi sources and theoretical estimates of the contribution of diffuse Galactic 
emission from the tot al sky i ntens ity. 

In a recent paper, Inoud ( 2011 ) estimated the contribution of unresolved 
7— ray-loud radio galaxies (GRLRG) to the GeV DGB. Their contribution 
is given by an integral similar to that used for calculating the EBL [see eq. 
(T5])], except that it requires the addition of an opacity term which takes the 
propagation of the 7— ray photons through the EBL into account. 

Calculating the contribution of GRLRG to the DGB requires knowledge of 
their intrinsic spectrum and 7— ray-luminosity function (LF). Using a sample 
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Figure 14: Limits on the optical depth at z = 0.1, 0.3, 0.5, and 1.0 as determined by the 
EBL models are compared to observationally determined optical depths for select blazars. 
The optical depth to 3C 66A is still uncertain. See §7 for more details. 



of 10 radio galaxies from the Fermi catalog, llnoud (120111 ) derived a correlation 
between their 5 GHz and 0.1-100 GeV luminosities. Using this correlation 
and the 5 GHz LF, he then derived the 7— ray-LF for these objects. 

The results show that GRLRG can account for about 25% of the intensity 
of the DGB above 100 MeV. Their cumulative spectrum is characterized by a 
power law to energies of ~ 10 GeV, after which it exhibits a decline resulting 
from EBL absorption. However, the observed Fermi-LAT spectrum does 
not exhibit such decline. The lack of this absorption signature in the DGB 
spectrum therefore suggests that nearby 7— ray sources or distant sources 
with hard 7— ray spectrum must be the major contributors to the DGB at 
energies above ~ 30 GeV. 

The derived GRLRG spectrum was derived using a sample of only 10 
galaxies to calculate their 7— ray LF. The model also contains some additional 
free parameters that characterize their intrinsic spectra and their overall 
contribution to the population of all radio galaxies. Clearly more data are 
required to resolve the origin of the DGB and detect any signature of EBL 
absorption in its spectrum. 
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8. Is the Blazar Spectrum Determined in the Source? 



In deriving upper limits on the EBL from TeV observations, it was tacitly 
assumed that blazar spectra are produced in the sources and attenuated 
en route to Earth. The detection of medium redshift blazars at z ~ 0.2 
with hard 7— ray spectra has revived alternate models for the origin of their 
spectra. Since hard photons from these redshifts are easily attenuated, one 
of these models proposes that the observed blazar spectra contains secondary 
photons that are produced close to the observer, and therefore more likely 
to survive. 

High resolution optical spectra of the blazars PKS 0447-439 and PMN J0630- 
24 obtained with the CTIO and NTT observatories show a clear absorp- 
tion line around 6280 A which, when attributed to the Mg II 2795.5 A and 
2802.7 A doublet place s them at redshifts z > 1.246 and z > 1.238, re- 
spectively (ILandtl . |2012| ). PKS 0447-439 was detected as a very-high energy 
7— ray so urce with H.E.S. S., showing an energy spectrum from 200 GeV - 
1.5 TeV flZech et all lioTlh . 

The optical depth of 1 TeV photons at such large redshift is about 20 (see 
Figures [12] and [13]) , so that any flux of ~ 1 TeV photons is attenuated by a 
factor of ~ 10~ 9 . This would suggest that the intrinsic 1 TeV luminosity of 
this blazar is about 10 10 larger than that of Mrk 421! Any attempts to solve 
this problem by lowering the 7 — 7 opacity of the intervening intergalactic 
medium will require the unrealistic reduction of the EBL intensity below the 
lower limits determined by the IGL (see Figure [7]). 

It is therefore unlikely that primary 7— rays could have reached the Earth 
from such distance, suggesting that most of the 7— rays from this object must 
be secondary photons created by the int eraction of cosmic-ray protons at 



relati vely close distances to the observer ( IKusenkol . l2012t lAharonian et al 



20121 ). The secondary photons are generated by the interactions of the pro- 



tons with the CMB and the EBL. For the secondary photons to be detected, 
they have to be produced within a distance A 7 , the mean-free-path of the sec- 
ondary 7— ra ys, of the Earth, an d be deflected into the viewing angle of the 
telescope (e.g lStanev et al.l . |2000| ). Their angular deviation from the primary 
source direction must therefore be smaller than the point spread function of 
the telescope (~ 0.1°), requiring the protons to travel on a "straight" trajec- 
tory until the interaction region. This requires their gyro-radius to be much 
larger than the source distance, which sets a range ~ 10~ 17 — 10~ 15 G on t he 



line-of-sight intensity of the intergalactic magnetic field ( lAharonian et al. 
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20121 ) . While this value is significantly smaller than commonly accep ted up- 



per limits for the intergalactic magnetic field (jKronbergl . Il994l 120101 ). weak 
fields are not ruled out. 

Hadronic particle acceleration is at the heart of the cascade model. This 
acceleration mechanism may hold the key to understanding blazar jets, and 
has far reaching consequences in high energy ast rophysics. Hadronic jet mod- 
els were brought forward by ( Mannheim . 19931) around the time of the first 
detection of TeV photons from a blazar ( iPunch et al.L 119921 ) . Protons at en- 
ergies of 10 18 — 10 19 eV are capable of reaching the threshold for photopion 
production. The subsequent pion-induced cascade inside the jet generates 
primary 7— ray photons which make a substantial contribution to the blazer 
spectrum. An additional contribution to the 7— ray spectrum comes from 
proton synchrotron radiation, and the synchrotron emission from the sec- 
ondary pions and muons. 

The viability of hadronic 7— ray emission models has been questioned, 
since they seemed to require exceedingly large magnetic fields to accelerate 
the protons and to confine them to the small emission region inferred from 
the observations of short duration TeV flares. However, hadronic 7— ray 
emission models have been put to test through TeV 7— ray observations of 
nearby blazars for which EBL absorption is negligible. The acceleration of 
UHE protons along relativistic jets requires large magnetic fields, typically 
several 10 G. The emission of gamma rays through a hadron induced cascade 
has to occur within a small emission region because of the observed variability 
of TeV flares. The fi rst big flare reporte d from Mrk 421 exhibited sub-hour 



scale flux variability (IGaidos et al.l . Il996l ). Assuming a Doppler factor of 10, 



the short flux variations set a limit of < 10 13 m, in the comovin g frame of the 



relati vistic jet, on the size of the gamma-ray emission region (IGaidos et al. 



19961 ). A magnetic field strength of at least ~ 50 G is required to prevent 
10 19 eV ions and protons from escaping the small emission region (jDermerl . 
20121 ). Such large magnetic fields are still consistent with the substantial, 



but plausible, power requirements for the magnetic field energy, which is less 
than 0.1% of the Eddington luminosity. Consequently, the observation of 
short duration TeV flares does not pose a major obstacle to the viability of 
hadronic models. 

Finally, hadronic jet emission models are also attractive for explaining 
the origin of ultra-high energy (UHE: E > 10 17 eV) cosm ic rays that are 



detected by air shower arrays such as AUGER and HiRes (iKampertl . 12012 



Sokolsky for the HiRes Collaboration! . 120101 ). 
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However, the need to resort to hadronic acceleration model to specif- 
ically explain the observed flux of PKS 0447-439 is probably premature. 
The redshift de t ermin ation to this blazar has recently been challenged by 
Fumagalli et al.l (120121 ) who independently obtained the PKS 0447-439 spec- 
trum using the 6.5 m Magellan Telescopes with a high S/N ratio of ~ 150 
in the 6270-6300 A region of interest. They point out the existence of an 
atmospheric telluric absorption line at the wavelength of the claimed Mg II 
absorption line, which they show ed was also present in the spectrum of two 
standard stars. iPita et al.l (120121 ) also called into question the redshift mea- 
surement of PKS 0447-439, providing clear evidence that the conclusions in 
Landtl fl2012h are incorrect. As a result, the redshift of PKS 0447-439 is still 
undetermined. 

Hadronic jet models can play an important role in explaining the origin 
of the DGB. The f eatureless DGB pow er spectrum derived from the Fermi- 
LAT observations ( lAbdo et al.l . |2010eJ) does not exhibit the EBL absorption 
signature around > 30 GeV, that is expected from the contribution of cosmo- 
logically distant 7— ray sources such as GRLRG, star burst, or active galaxies 
(for a discussion see llnoud . l2011t llnoue &: Iokal . 120121 and references therein) . 
Hadron induced pair cascades interacting with the CMB can generate sec- 
ondary 7— rays that will con tribute to the GeV DG B, filling in the "missing" 
flux at energies > 30 GeV (lAharonian et al.l . 120121 ) . Further measurements 
extending the DGB to ~ 1 TeV will provide important constraints on the 
contribution of secondary 7— rays from blazars. 



9. EBL Limits on "Exotic" Energy Releases in the Universe 

The intensity and spectral shape of the EBL contain the memory of 
all energy releases in the universe since the epoch of recombination. The 
COBE/DIKBE limits on the UV to near-IR regions of the EBL proved use- 
ful in ruling out various "exotic" sources of energy in the early u niverse, such 



as de caying particles, exploding stars, or very massive objects ( jDwek et al. 



19981 . and references therein). 



Primordial (Population III) stars were suggested by ISalvaterra fc Ferrara 



( 120031 ) as the sourc e of the excess 1-5 / im di ffuse emissio n above the IGL in - 
tensity detected by lMatsumoto et al.l ( 120051 ). However, ( jDwek et al.l . l2005bl ) 
ruled out an extragalactic origin for the excess emission, since it would have 
produced a physically unreal istic intrinsi c 7— ray spectrum of the blazar 
PKS 2155-304. Furthermore, iDwek et al.l (l2005al ) showed that such origin 
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Figure 15: The contribution of dark stars to the EBL for two different Limits on the EBL 
intensity. Lower limits, blue squares) are determined by the intensity of the IGL. Upper 
limits (red circles) are determined by absolute measurements of the EBL. The data used 
in the figure are listed in Tables 3-5 in bold. The shaded area depicts the range of the 
allowed EBL intensity as determines by UV to sub millimeter observations. 



would have required a Pop III star formation rate and energy output to be 
significantly higher than that pre dicted by hierarchical m odels for structure 
formation in a ACDM universe ( Bromm &: Loebl. 2002 ). Using theoretical 
limits on their formation rate, ( Dwek et al. . 2005al ) concluded that Pop III 
stars can contribute only a fraction of the EBL intensity. 

A more detailed study o f the contri b ution of Pop III stars to the EBL 
intensity was conducted by iRaue et al.l (120091 ) . Using 7— ray derived con- 
straints on the EBL intensity they set a limit of 0.3 to 3 M Q Mpc -3 yr" 1 on 
the formation rate of Pop HI star s in th e z = 7 — 14 redshift interval. 



More recently, ( IMaurer et al 



20121 ) used EBL limits to constrain the 
properties of dark stars (DS) in the early universe. Dark stars are objects 
that have either accreted or captured weakly interacting massive particles 
(WIMPs, a dark matter candidate), which by annihilating inject energy into 
the stars before their radiative output is dominated by standard nuclear 
fusion processes. The formation rate of these stars, their luminosity and 
spectrum, and their effective lifetime are all free parameters of the model. 
Figure [TBI depicts the EBL, to which the contribution of dark stars was added 
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to the IGL intensity, for two different sets of parameters. The blue curve rep- 
resents the EBL with a contribution of 106 M Q DS with a surface temperature 
of 5,000 K, and the red curve that with the added contribution of 690 M Q DS 
with a surface temperature of 7,500 K. The colder dark stars are obviously 
ruled out, but the hotter ones are marginally consistent with current EBL 
limits, dominating the intensity of the IGL in the ~ 2 — 10 fim wavelength 
region. 

10. Summary and Future Direction 

Very high energy 7— rays emitted from extragalactic sources are attenu- 
ated en route to earth by 7 — 7 interaction with EBL photons. 7— ray obser- 
vations can therefore be used to set limits on the EBL intensity, provided that 
the intrinsic 7— ray spectrum of the sources is known. Conversely, knowledge 
of the EBL can be used to determine the intrinsic 7— ray spectrum of the 
different sources, thereby provide important constraints on mechanisms for 
their production. The main issues and results discussed in this review can 
be briefly summarized as follows: 

1. The EBL spectrum consists of two broad peaks, one at A ~ 1 /im, 
representing the cumulative gravitational and nuclear energy releases 
by stars and AGNs over cosmic history. Their energy has been partially 
absorbed by dust and reradiated at IR wavelengths. This thermal 
dust emission component generates a second peak at ~ 100 — 200 /xm. 
Current limits and detections of the EBL were presented in Figure [7J 
and the partitioning of its total intensity into the different emission 
components was presented in Table 6; 

2. The cross section for the 7—7 interaction is broad and peaks at energies 
E 7 (TeV) ~ 0.86A(/xm). Consequently, the energy dependence of the 
7— ray opacity, a product of the cross-section with the number density 
of EBL photons, reflects the spectral variation in the EBL. We identified 
three potential breaks in the spectrum of 7— ray sources, one occurring 
between 10 and 500 GeV, a second break at 1 TeV and third at 10 TeV 
(see Figures [3] and H]) ; 

3. The spectral breaks at GeV and TeV energies have been used to set 
limits on the EBL intensity at near-IR wavelengths, and on the relative 
intensities of the peak of the stellar emission and the trough between 
the stellar and dust emission components of the EBL; 
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4. 7 — ra y derived limits on the EBL vary as different studies used different 
assumptions on the intrinsic 7— ray source spectra (see Figure ITU]) . The 
strictest limits on the EBL are around the 1 and ~ 10 fim wavelength 
regions of the EBL, with some approaching the lower limit on the EBL 
intensity set by the IGL (see Figure [XT]) ; 

5. Recent IR space and ground-based observations have resulted in closer 
agreement between the EBL limits and detections derived from mea- 
surements of the absolute sky brightness and lower limits set by the 
integrated light from galaxies. A summary of the recent observational 
status of the EBL is presented in Tables 3-5. The major gap in our 
knowledge of the EBL is in the A ~ 10—70 //m wavelength region, where 
the thermal emission from interplanetary dust (the zodiacal light) dom- 
inates the brightness of the sky. This wavelength region is best probed 
by nearby radio and starburst galaxies at energies > 10 TeV; 

6. Models of the EBL, employing various methods for determining the 
evolution of the galaxies' spectral energy density with redshift seem 
to agree on the intensity and spectral shape of the UV to near-IR 
component of the EBL. However, there are considerable differences in 
their treatment of the redistribution of the intrinsic stellar and AGN 
output at IR wavelengths (see Figure |9]); 

7. The discovery of TeV blazars with alleged redshifts z > 1 has revived 
an alternative model for the creation of these high energy photons. In 
this model hadronic jets produce a cascade of ~ TeV 7— ray photons en 
route to earth, circumventing the attenuation problem from such high 
redshift sources (see §8). However, the redshift of these sources has 
been disputed, so the need to resort to such hadronic model is highly 
premature; 

8. Dark matter has been invoked to postulate the existence of a new class 
of primordial stars, powered by dark matter annihilation, instead of 
nuclear fusion. The energy release from these so-called dark stars could, 
in principle, lead to an observable signature in the EBL spectrum (see 
Figure [T5]) ; and finally 

9. Observations suggest that the universe is essentially transparent to 
7— rays with energies < 2 TeV up to z ~ 0.2, and energies of < 400 GeV 
up to z ~ 0.4 (see Figure [T4"j) . 

The future prospects of intensified EBL studies with 7— rays are promis- 
ing, especially when considering the progress made over the last 5 years 
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through the operation of Fermi and the current generation of atmospheric 
Cherenkov telescopes. These instrument yielded combined energy spectra for 
at least 3 dozen extragalactic sources that cover up to 5 orders of magnitude 
in energy. These data have provided a first glimpse of the transition region 
from a transparent universe at 1 GeV, to TeV energies at which the universe 
gradually turns opaque with increasing redshift. 



Furthermore, Fermi detected a large set of blazars (lAckermann et al. 



201 ll ) including 310 FSRQs, 395 BL Lacertae objects and 156 candidate 



blazars, raising the specter for EBL studies with sizable source samples and 
different sources classes with CTA. Large samples of blazars have the po- 
tential to constrain the EBL in the optical/near-IR and mid- IR through a 
better understanding of the blazar subclasses (FSRQs, LBL, IBL, HBL) and 
their intrinsic spectra, and better photon statistics for the measurement of 
the redshift dependence of any spectral feature attributable to the EBL. 

Recent discoveries of TeV emission from nearby radio- and starburst 
galaxies also gives rise to future prospects of extending the reach of 7— ray 
constraints up to the far-IR through energy spectra spanning a few GeV up 
to 100 TeV. These will be important to extend measurements of the 7— ray 
opacity imposed by the EBL all the way from UV/optical/near-IR/mid-IR 
to the far-IR and thus provide additional constraints to the intensity ratios 
between the different wavelength regimes of the EBL. 

Limitations of precision EBL studies with 7— rays arise from technical 
reasons, but are not unsurmountable. First, systematic uncertainties in the 
measurements of 7— ray spectra with atmospheric Cherenkov telescopes arise 
from uncertainties in modeling the Earth's atmosphere, which translates into 
uncertainties in the attenuation of Cherenkov light from electromagnetic cas- 
cades. To first order this affects the absolute energy scale of the measured 
7— ray energies and is typically quoted at a level of 15% to 20%. Effects 
on the spectral index and shape are generally of second order but require 
detailed studies. Additional uncertainties lie in the absolute instrument cal- 
ibration of Cherenkov telescopes, e.g., mirror reflectivity, light losses in the 
focal plane, and uncertainties in the quantum efficiency of the photodetec- 
tors, affecting the light throughput and absolute energy scale. Uncertainties 
in the spectral indices and spectral shape are typically quoted at the level of 
Ar = 0.05 — 0.1 or better, and depend mostly on the 7— ray selection effi- 
ciency derived from Monte Carlo simulations and detailed detector modeling. 
Currently most energy spectra published (except for few flaring sources) are 
dominated by statistical uncertainties and leave much room for improvement 
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through a more sensitive instrument such as CTA with a collection area of 
~ km 2 for 100 GeV to a few TeV energies, and several km 2 for higher ener- 
gies. At the same time calibration techniques are continually improving to 
provide a better handle on systematics as well. 

Astrophysical challenges facing EBL studies with 7— rays arise from the 
fact that currently o nly 55% of all Fermi det ected BL Lacs have reliable red- 
shift measurements ( Ackermann et al. . 201 ll ). Attempts to increase the red- 
shift identifications through dedicated optical follow-up observations of the 
host galaxies of Fermi detected blazars during low flaring states are promis- 
ing to increase the fraction of blazars with known redshift. The redshift 
distribution for FSRQs detected by Fermi peaks at z = 1 extending to z = 
3.1, while BL Lacs peak at z = 0.2 reaching up to a redshift of z = 1.5. 
As already indicated, GeV/TeV detections of FSRQs are promising to allow 
constraints to the EBL in the UV/optical. Furthermore, the first detection 
of FSRQs by atmospheric Cherenkov telescopes out to redshifts of z = 0.5 
provides the potential to complement EBL studies with a different class of 
sources, with different intrinsic spectral properties than BL Lacs and, most 
importantly, whose redshifts are known. 

Detection and limits on the EBL have been obtained by direct absolute 
measurements, and by galaxy number counts that provided lower limits on 
its intensity. Direct measurements are aided by ground-based or space-based 
observations that resolve the foreground emission from Galactic stars. Fu- 
ture observations can considerably improve direct absolute measurement of 
the EBL by determining the absolute brightness of the zodiacal light from 
high resolution observations of reflected solar Fraunhofer lines, or by making 
absolute sky measurements from the outer solar system, thus effectively re- 
moving any foreground contribution from the zodiacal cloud. Lower limits on 
the EBL can converge to the EBL itself with observations of sufficient depth 
and resolution to resolve all the galaxies that contribute to its intensity. 

The James Webb Space Telescope (JWST) is a large (6.6 m) IR space ob- 
servatory, passively cooled to temperatures below 50 K, that will be launched 
into orbit at the second Earth-Sun Lagrange point (L2). The wide band ob- 
servations with NIRCam at 2 - 5 /im and MIRI at 6 - 25 /im will fill in crucial 
gaps in our knowledge of the EBL intensity at these wavelengths. NIRCam 
will cover a 10 arcmin 2 field of view (FOV) and will resolve galaxies down 
to a 10<7 flux limit of 11 nJy at 2 /j,m in 10,000 s, making it about 10 3 more 
sensitive than the Spitzer/IRA C 3.6 um filter . Wit h unparalleled resolu- 
tion and sensitivity, the JWST ([Gardner et al.l . 120061 ) will resolve the EBL 
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at near-IR wavelengths. The MIRI instrument operating in its broadband 
imaging mode at 10 and 21 /iin will have a smaller FOV of ~2.6 arcmin 2 , 
and sensitivities of 700 nJy and 9 /zJy, respectively. It will resolve galaxies 
down to the confusion limit which, because of its large telescope size, will be 
significantly fainter than that of the Spitzer/MIPS instrument at 24 /im. 

Absolute measurements of solar absorption line profiles in the zodiacal 
light, for example the Fraunhofer Mg I line at 1.182819 /zm or the Si I line 
at 1.210354 fim can be used to determine and remove the contribution of 



the ZL to t he total sky brightness in the J (1.25 fim) band (IKutyrev et al. 



2004] . [2008|). High resolution (R ~ 20,000) observations of these lines us- 



ing the ground-based spectrometer ZEFIR (Zodiacal Emission determination 
through Fraunhofer IR lines) will determine their equivalent width (EW) in 
the ZL spectrum. The continuum contribution of the ZL in the J-band can 
then be determined from knowledge of the EW of these lines in the solar 
spectrum. The DIRBE instrument provided absolute measurements of the 
sky brightness in this band, and ground-based sky surveys (2MASS) have re- 
moved most of the contribution of Galactic stellar emission from this band. 
Measurements of the absolute brightness of the ZL will thus provide the 
absolute intensity of the EBL in the DIRBE 1.25 fim waveband. Similar 
measurement could also be made from the ground in the 2.2 /zm band. 

Finally, the contribution of the ZL to the observed sky brightness can be 
largely eliminated by mapping the absolute intensity of the diffuse emission 
from the sky from beyond the zodiacal cloud. A satellite mission conducting 
such EBL measurements from the outer solar system will be an important 
step in that direction (Cooray, 2011). 

The next decade will see significant improvement in 7-ray technology, and 
measurements of the EBL, providing new insights into the spectrum of the 
different 7-ray sources and the mechanism that generate them, and into the 
spectrum of the EBL which will provide new constraints on the history of 
all nuclear and gravitational energy releases in the universe. The synergy 
between TeV 7-ray astronomy and EBL research will ensure that any devel- 
opment in each field will greatly benefit both. 
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Table 1: Glossary of Abbreviations of Spacecrafts 1 , Telescopes, and Instruments 


Abbreviation 


Full Name 


Akari 


Infrared imaging satellite (ASTRO-F) 


BLAST 


Balloon-borne Large-Aperture Submillimeter Telescope 


COBE 


Cosmic Background Explorer 


DIRBE 


Diffuse Infrared Background Experiment 


FIRAS 


Far Infrared Absolute Photometer 


CTIO 


Cerro Tololo Inter-American Observatory 


GALEX 


Galaxy Evolution Explorer 


Herschel 


Herschel Space Observatory 


PACS 


Photodetector Array Camera 


SPIRE 


Spectral and Photometric Imaging Receiver 


HST 


Hubble Space Telescope 


WFPC2 


Wide Field Planetary Camera 


NICMOS 


Near IR Camera and Multi-Object Spectrometer 


IRTS 


Infrared Telescope in Space 


ISO 


Infrared Space Observatory 


ISOCAM 


ISO Camera 


JCMT 


James Clerk Maxwell Telescope 


SCUBA 


Submillimeter Common User Bolometer Array 


NTT 


New Technology Telescope 


Spitzer 


Spitzer Space Telescope 


IRAC 


Infrared Array Camera 


MIPS 


Multiband Imaging Photometer 


Subaru 


Optical, near-IR telescope 


Pioneer 


Interplanetary spacecraft 


WMAP 


Wilkinson Microwave Anisotropy Probe 


2MASS 


Two Micron All Sky Survey 


CTA 


Cherenkov Telescope Array 


Fermi 


Fermi Gamma-Ray Space Telescope 


H.E.S.S. 


High Energy Stereoscopic System 


IACT 


Imaging Air Cherenkov Telescope 


MAGIC 


Major Atmospheric Gamma-Ray Imaging Cherenkov Telescope 


Milagro 


Gamma-ray and cosmic-ray telescope 


VERITAS 


Very Energetic Radiation Imaging Telescope Array 



1 Spacecraft names are presented in italics, and their instruments are in- 
dented. 



Tabic 2: Extragalactic 7— ray sources with GeV and TeV spectral information. 



Name 


Class 


rcdshift 


^GeV 




Range 


[TeV] 


References 


Centaurus A 


Radio 


0.0008 


2.76+0.05 


2.7+0.5 


0.2 


- 5 


[1], M 


M82 


SB 


0.00085 


2.2+0.2 


2.5+0.6 


0.7 


- 4 


[3], [4] 


NGC253 


SB 


0.00093 


1.95+0.4 


2.24+0.14 


0.3 - 


50 


[3], [5] 


M87 


Radio 


0.0036 


2.17+0.07 


2.5+0.2 


0.2 - 


10 


[6], [7], 8], [9] 


NGC 1275 


Radio 


0.018 


2.00+0.02 


3.96+0.37 


0.1 - 


0.3 


[10], [11] 


IC 310 


Radio 


0.0188 


2.10+0.19 


2.0+0.14 


0.1 


- 7 


[12], [13], [14] 


Markarian 421 


HBL 


0.031 


1.77+0.01 


2.48+0.03* 


0.1 


- 5 


[15] 


Markarian 501 


HBL 


0.034 


1.74+0.03 


2.51±0.05 A 


0.1 - 


10 


[16] 


1ES 2344+514 


HBL 


0.044 


1.72+0.08 


2.78±0.09 A 


0.3 


- 2 


[6], [17] 


Markarian 180 


HBL 


0.046 


1.74+0.08 


3.3+0.70 


0.2 


- 1 


[6] [18] 


1ES 1959+650 


HBL 


0.047 


1.94+0.03 


2.72+0.14 


0.2 


- 2 


[6], [19] 


AP Lib* 


LBL 


0.048 


2.05+0.04 


2.5+0.2 


0.3 


- 2 


[6], [20] 


BL Lacertae 


LBL 


0.069 


2.11+0.04 


3.6+0.5 


0.2 


- 1 


[6], [21] 


PKS 2005-489 


HBL 


0.071 


1.78+0.05 


4.0+0.4 


0.2 


- 2 


[6], [22] 


W Comae 


IBL 


0.103 


2.02+0.03 


3.81+0.35 


0.3 


- 1 


[6], [23] 


PKS 2155-304 


HBL 


0.116 


1.84+0.02 


3.53+0.05 


0.4 


- 5 


[6], [24] 


B3 2247+381 


HBL 


0.119 


1.84+0.11 


3.2+0.5 


0.2 


- 1 


[6], [25] 


RGB J0710+591 


HBL 


0.125 


1.53+0.12 


2.69+0.26 


0.3 - 


4.6 


[6], [26] 


H 1426+428 


HBL 


0.129 


1.32+0.12 


3.50+0.35 


0.3 - 


10 


[6], [27] 


1ES 1215+303 


IBL 


0.13^ 


2.02+0.02 


2.99+0.15 


0.1 


- 1 


[6], [28] 


1ES 0806+524 


HBL 


0.137 


1.94+0.06 


3.6+1.0 


0.3 - 


0.7 


[6], [29] 


1RXS J101015. 9-311909 


HBL 


0.143 


2.24+0.14 


3.14+0.53 


0.3 


- 1 


[6], [20] 


1ES 1440+122 


IBL 


0.163 


1.41+0.18 


3.3+0.7 


0.3 


- 1 


[6], [30] 


H 2356-309 


HBL 


0.165 


1.89+0.17 


3.09+0.24 


0.3 


- 2 


[6], [31] 


VER J0648+152 


HBL 


0.179 


1.74+0.11 


4.4+0.8 


0.3 - 


0.8 


[6], [32] 


1ES 1218+304 


HBL 


0.184 


1.71+0.07 


3.07+0.09 


0.2 


- 2 


[6], [33] 


1ES 1101-232 


HBL 


0.186 


1.80+0.21 


2.88+0.17 


0.16 


- 3.3 


[6], [31] 


RBS 0413 


HBL 


0.19 


1.55+0.11 


3.18+0.68 


0.25 


- 1 


[6], [33] 


PKS-0447-439 


HBL 


0.205 


1.86+0.02 


4.36+0.49 


0.25 


- 1 


[6], [34] 


1ES 1011+496 


HBL 


0.212 


1.72+0.04 


4.0+0.50 


0.25 


- 0.6 


[6], [35] 


1ES 0414+009 


HBL 


0.287 


1.98+0.16 


3.44+0.27 


0.25 


- 1.2 


[6], [36] 


S5 0716+714 


LBL 


0.31 


2.01+0.02 


3.45+0.54 


0.25 


- 1.2 


[6], [37] 


1ES 0502+675 


HBL 


0.416* 


1.49+0.07 


3.92+0.35 


0.25 


- 1 


[6], [38] 


4C 21.35 


FSRQ 


0.43 


2.12+0.02 


3.75+0.27 


0.07 


- 0.4 


[6], [39] 


3C 66A 


IBL 


0.44* 


1.85+0.02 


4.1+0.4 


0.22 - 


0.45 


[6], [40] 


3C 279 


FSRQ 


0.536 


2.22+0.02 


3.03+0.9 


0.1 - 


0.35 


[6], [41] 



* Spectrum is not well fit by a powcrlaw. A Spectrum shows variations. £ redshift uncertain. ^? redshift 
recently given by Abdo et al. arXiv:1108.1420vl, different redshift was considered viable by Colin et al. 
2011 

[l jAbdo et al] feOlOdl): [21 [Aharonian et alj (|2009T ): [3] I Abdo et alj j2010d1; [41 lAcciari et al] 
for t he H.E .S.S. collaboration); [6 1 lAckermann et al] J201lf>; 
lAcciari et all d2010bl1: [91 I Aharonian et at] ll2006bh; [1 ] I Abdo et al] 



llAbramowski et a l | |201:' 



Referen ces: 
J2009bft: [51 . 

171 iBerger et alj feoilaf): [81 

j2009t); [111 iHildebrand et alj (201 ifl: [1 21 iNeronov et al] ll2010l) ; [13] lAleksic et al] ifeOlObT) : [14 1 
lAckermann et alj J2011T); [151 lAbdo et alj [feoilcT); [16 1 lAbdo et al] ll2011aft;~il7l lAcciari et al] d2011bt>: 
[181 lAlbert et alj d2006al1: [191 lAlbert et alj (l20d6bh : [20lCerruti M. for the H.E.S.S. Collaboration I 
d201lT); [211 lAlbert et al] <2007ah: [221 lAharoriian et al] <2005ah; [231 lAcciari et al] J200ST ) ; [24 ] 
lAharonian et all ll2007bl 1: [ 25flBerger et a I teOllbh : [ 261 lAcciari et alj J2010al1: [271 iPetrv etldTll2002ft : 
[28llColin et aljfeOllD; [2 91 lAcciari et alj d2009al'): ~[3 0] iBenbowl ||2011|1; [3 l TlAharonian et al] d2006aT) ; [32 1 
lErrando. M. et al.l j201lT) [331 lAcciari et all d2010j> : [3 31 ISenttirk et alj llpbllh : [341 IZech et alj j201lfi: 
[351 lAlbert et al] j2007bh ; [361 IVolpe et al] 1120111); [371 lAnderhub et alj j2008h; [381 [Benbow, W. et al] 
d201ir> ; [39] lAleksic et all fcOllaT ); [401 lAcciari etall d2009cf ); [411 lAleksic : etldl j2011bil 
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Table 3: Limits and detection of the extragalactic background light (EBL). 
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3.6 






5piteer/IRAC 


[16] 


8.0 


2.6 






S'prfzer/IRAC 


[16] 



ABS=absolutc measurement; IGL=intcgratcd galactic light; STK=lower limits from stacking analysis; 
EXT=extrapolated intensity from dN/dS 

1 Calculated for a 2.2 fim inte nsity of 20.0 nW m ~ 2 sr" 1 an d the Kelsall et al. (1 998) ZL model 
Ref erences: [11 IXu et all J2005D; [21 iGardner et al] (|2000h ; [3] iMilliard et ahl jl992h; [41 iBernsteid j2007tl ; 
[51 iTotani et all j200ll) ; [6 1 iMadau fc Pozzettil I2OO0I); [71 iMattilal <199CMLeinert et all dl998fl ; [8] 
Matsuoka et all ir201lT) : [9] ll)ube et all dl979TllLe"inert et all lll998l1: [ 101 iLevenson et all (1200711: [111 
Cambresv et all j200ll1; [12~TKeenan et all ll201Cft; [131 iMatsumotd ll200ll); [141 iLevenson fc Wright! 
j2008fl ; [15l lDwek fc Arendtj Jl998f ): [16l lFazio et all ll2004h ; [171 lArendt fc DweU j20o7J) 
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Table 4: Limits and detection of the extragalactic background light (EBL). 
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EXT 
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Comment 


Reference 


15 


2.4 ±0.5 
2.7±0.6 
1.9 ±0.5 








KO/ISOCAM 
JSO/ISOCAM 
Akari 


Elbaz ct al. (2002) 
Metcalfe et al. (20031 
Hopwood ct al. (2010) 


16 


2.2 ±0.2 
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Tcplitz ct al. (2011) 


24 


1.9 ±0.6 
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Papovich ct al. (2004) 




1.8 ±0.2 
2.29 ±0.09 
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Charv et al. ('2004) 
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Matsuura ct al. (2011) 
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Akari 


Matsuura ct al. (2011) 
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8.35 ±0.95 


9.4± 1.1 




< 34 (22 ±6) 
12.5 ± 5 

24.6 ±2.5 ±8 
23.4 ±6.3 
14.4±6.0 

6.6 ±1.8 ±2.1 


COBE/DIRBE (D, KZL) 
COBE/DIKBE (D, WZL) 

COBE/DIKBE 
COBE/DIKBE (D, KZL) 
DIRBE (F, WZL) 
Herschel/PACS 
Spitzer/MlPS 


Hauscr ct al. (1998) 
Wright (2004) 
Finkbeiner et al. (2000) 
Laeache et al. (2000) 
Dole et al. (2006) 
Berta et al. (2010) 
Pcnin ct al. (2012) 


140 








25.0 ±6.9 
15.0 ±5.9 

32 ± 13 
24.2 ± 11.6 
22 ±7 
20.1 ±3.4 ±1.1 
12.4 ±6.9 
12.6±6.0 


COBE/DIKBE (D, KZL) 
COBE/DIRBE (F, KZL) 

COBE/DIKBE 
COBE/DIKBE (D, KZL) 
COBE/DIKBE (D, WZL) 
Akari 

COBE/DIKBE (F, WZL) 
COBE/EIRAS 


Hauscr ct al. (1998) 
Odegard et al. (2007) 
Schleeel et al. (1998) 
Laeache et al. (2000) 

Wright (2004) 
Matsuura et al. (2011) 
Dole et al. (2006) 
Fixsen et al. (1998) 
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Table 5: Limits and detection of the extragalactic background light (EBL). 
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A(//m) 


IGL 


STK 


EXT 


ABS 


Comment 


Reference 


i fin 


9.49 ± 0.59 
8.9 ± 1.1 


11.4 ±0.7 


14 6+ 7 - 1 


13.7 ±3.9 ±0.8 

14.4 ±0.8 ±2.3 
13.7±6. 1 


j7 eVSCflcL/ rfiba 

Hersc/ieZ/PACS 
Spiteer/MIPS 
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uoie cz ai. \ zuud ) 
Matsuura et al. (2011) 
Berta et al. (2010) 
Bethermin et al. (2010a) 
Penin et al. (2012) 
Fixsen et al. (1998) 
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19.1 ±5.6 ±5.3 








ISOPHOT 


Juvcla ct al. (2009) 


240 








13.6 ±2.5 
13 ± 2.5 

12.7 ± 1.6 
17 ±4 

11.0 ±6.9 
12.3 ± 2.5 
10.9±^.3 


CO-RE/DIRBE (D, KZL) 
CO.RE/DIRBE (D, WZL) 
CORE/DIRBE (F, KZL) 

CORE/DIRBE 
COBS/DIRBE (D, KZL) 
CORE/DIRBE (F, WZL) 

COBE/FIRAS 


Hauscr ct at (1998) 

Wright (2004) 
Odeeard et al. (2007) 
Schleael et al. (1998) 
Laeache et al. (2000) 
Dole et al. (2006) 
Fixsen ct al. ( 1998) 


250 


1.73 ±0.33 

ice in o/"i 

1.55 ± 0.30 


8.6 ±0.6 
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7.40 ± 1.42 
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II, _ 11 / CTDTTDTI 1 
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Oliver et al. (2010) 
Bcthcrmm ct ai. (2012 J 

Fixsen et al. (1998) 


350 


06+ 05 
u ' UD -0.04 

0.63 ±0.18 
0.77 ±0.16 


4.93 ± 0.34 
9 S+ 1 - 8 

4.50 ±0.90 




3.6±2.1 


BLAST 

BLAST 
Hersc/ieZ/SPIRE 
tfersc/ieZ/SPIRE 

FIRAS 


Marsden et al. (2009) 
Bethermin et al. (2010b) 

Oliver et al. (2010) 
Bethermin et al. (2012) 

Fixsen ct al. (1998) 


500 


m+ 01 

uul -0.01 

0.15 ±0.07 
0.14 ±0.03 


2.27 ±0.20 
1 4+ 2 - 1 

1.54 ±0.34 




2.4±0.9 


BLAST 

BLAST 
Herschel/SPIRE 
Herschel/SPIRE 

FIRAS 


Marsden et al. (2009) 
Bethermin et al. (2010b) 

Oliver et al. (2010) 
Bethermin et al. (2012) 

Fixsen ct al. (1998) 


850 


0.12 ±0.03 
0.24 ±0.03 






0.5±0.21 


SCUBA 
SCUBA 
COBE/FTRAS 


Coppin ct al. (2006) 
Zemcov et al. (20101 
Fixsen ct al. (1998) 


200-1000 








a(±f) k vB„{T) 


FIRAS 


Fixsen et al. (1998) 1 


1 a = (1.3± 


= 0.4) x 1CT 5 ; 


fc = 0.64 ±0.12; T = 


(18.5 ± 1.2) K- A 


= 100 A*m 
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Table 6: Total EBL Intensity of Different Models (nW tcT 2 sr" 1 ) 



Model 


Stars 1 


Dust 2 


Total 


Franceschini et al. (2008) 


25 


40 


65 


Dommeuez et al. (2011) 


25 


44 


69 


Gilmore et al. (2011) 


25 


23 


48 


Mazin & Raue (2007) 


30 


26 


56 


Finke et al. (2010) 


27 


20 


47 


Stecker et al. (2006) 


61 


35 


96 


Observations 3 


23-93 


20-110 


42-202 


From CSFR 4 






21-66 



1 Integrated intensity from 0.1 to 10 jum. 

2 Integrated in tensity from 10 to 1000 /im. 

3 From limits and detections of the EBL (Figure [7]) 

4 Total integrated intensity inferred from the CSFR (FigjSJ) 
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